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ABSTRACT 

We derive the basic properties of seven Galactic open clusters containing Cepheids and construct their period- 
luminosity (P-L) relations. For our cluster main-sequence fitting we extend previous Hyades-based empirical 
color-temperature corrections to hotter stars using the Pleiades as a template. We use BVIcJHK s data to test 
the reddening law and include metallicity effects to perform a more comprehensive study for our clusters 
than prior efforts. The ratio of total to selective extinction Ry that we derive is consistent with expectations. 
Assuming the LMC P-L slopes, we find (My) =-3.93 ±0.07 (statistical) ±0.14 (systematic) for 10-day period 
Cepheids, which is generally fainter than those in previous studies. Our results are consistent with recent HST 
and Hipparcos parallax studies when using the Wesenheit magnitudes W(VI). Uncertainties in reddening and 
metallicity are the major remaining sources of error in the V-band P-L relation, but a higher precision could 
be obtained with deeper optical and near-infrared cluster photometry. We derive distances to NGC 4258, the 
LMC, and M33 of (m-M) = 29.28 ±0.10, 18.34 ±0.06, and 24.55 ±0.28, respectively, with an additional 
systematic error of 0.16 mag in the P-L relations. The distance to NGC 4258 is in good agreement with the 
geometric distance derived from water masers [A(m— M)o = 0.01 ± 0.24], our value for M33 is less consistent 
with the distance from an eclipsing binary [A(m-M)o = 0.37 ±0.34], and our LMC distance is moderately 
shorter than the adopted distance in the HST Key Project, which formally implies an increase in the Hubble 
constant of 7% ±8%. 

Subject headings: Cepheids — distance scale — galaxies: individual (M33, NGC 4258) — Magellanic Clouds 
— open clusters and associations: general — stars: distances 



1. INTRODUCTION 

The Cepheid period-luminosity (P-L) relation has played a 
key role in the determination of distances within the Local 
Group and to nearby galaxies. By extension, it is also crucial 
for the calibration of secondary distance indicators used to 
determine the Hubble constant: see, for example, the Hubble 
Space Telescope (HST) Key Project on the extragalactic dis- 
tance scale (Freedman et al. 2001, hereafter HST Key Project) 
and the HST program for the luminosity calibration of Type 
la supernovae (SNe la) by means of Cepheids (Sandage et al. 

2006) . The cosmic distance scale was usually established by 
defining the P-L relations for Cepheids in the Large Magel- 
lanic Cloud (LMC) because of its numerous Cepheids, many 
of which have been discovered as a result of microlensing 
campaigns (e.g., Udalski et al. 1999a). Despite decades of 
effort, however, there have been persistent differences in the 
inferred LMC distance from different methods (see Benedict 
et al. 2002a), including ones involving the same basic calibra- 
tors such as RR Lyrae or Cepheid variables. 

The Galactic P-L relationship, on the other hand, was tra- 
ditionally established using open clusters and associations 
containing Cepheids (see Feast & Walker 1987; Feast 2003; 
Sandage & Tammann 2006). There is also a steadily grow- 
ing body of parallax work from the Hipparcos mission (Feast 
& Catchpole 1997; Madore & Freedman 1998; Lanoix et al. 
1999; Groenewegen & Oudmaijer 2000; van Leeuwen et al. 

2007) and the HST (Benedict et al. 2002b, 2007), as well 
as distances derived from the Baade-Becker-Wesselink mov- 
ing atmospheres method (Gieren et al. 1997, 1998, 2005) and 
interferometric angular diameter measurements (Lane et al. 
2000, 2002; Kervella et al. 2004a,b). In the current paper we 
focus on the Cepheid distance scale as inferred from Galactic 



open clusters and its applications for the extragalactic distance 
scale. 

There are strengths and weaknesses in all of the methods 
used to establish the Galactic Cepheid distance scale. The 
absolute calibration of the P-L relation requires not only ac- 
curate distance measurements but also appropriate account- 
ing for the effects of interstellar extinction and reddening be- 
cause Galactic Cepheids are heavily obscured with an average 
E(B-V) of order 0.5 mag (Fernie 1990). Reddening can be 
inferred more precisely for clusters, while field Cepheids are 
more numerous. This difficulty can be partially overcome by 
the usage of the "reddening-free" or Wesenheit index (Freed- 
man et al. 1991), but there are embedded assumptions about 
the extinction law even in such a system. 

We believe that the time is right for a systematic reappraisal 
of the Cepheid distance scale as inferred from open clusters. 
Recent parallax work employing the Wesenheit index (Bene- 
dict et al. 2007; van Leeuwen et al. 2007) has yielded con- 
sistent claims of a somewhat smaller LMC distance modulus 
(« 18.40 mag) than the ones adopted by the HST Key Project 
(18.50 mag) or Sandage et al. program (18.54 mag). In addi- 
tion, there are now independent geometric tests of distances to 
other galaxies, which have large numbers of Cepheids. These 
include massive eclipsing binaries in systems such as M31 
(Ribas et al. 2005) and M33 (Bonanos et al. 2006) and astro- 
physical water masers in NGC 4258 (Herrnstein et al. 1999). 
As described below, the availability of both more complete 
data in various photometric passbands and better theoretical 
templates permits a more comprehensive look at the open 
cluster Cepheid distance scale, which can in turn be com- 
pared with the alternate methods described above. We can 
also use the open clusters with Cepheids to test our ability to 
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derive precise distances and reddenings to heavily obscured 
and poorly studied systems. We believe that the net result is 
a more secure and robust determination of the extragalactic 
distance scale. 

We have been engaged in a long-term effort to create stellar 
evolution models with the latest input physics and to generate 
isochrones that are calibrated against photometry in local star 
clusters with accurate distances. After verifying that the mod- 
els are in agreement with the physical parameters of binaries 
and single stars in the Hyades (Pinsonneault et al. 2003, 2004, 
hereafter Papers I and II, respectively), we developed a pro- 
cedure for empirically correcting the color-effective tempera- 
ture (T e ff) relations to match photometry in the Hyades (Paper 
II). Isochrones generated in this way accurately reproduce the 
shapes of the main-sequence (MS) in several colors, allowing 
the determination of distances to nearby clusters to an accu- 
racy of <~ 0.04 mag in distance modulus (or 2% in distance) 
(An et al. 2007, hereafter Paper III). 

We extend the Hyades-based empirical corrections on the- 
oretical isochrones using the Pleiades and apply our tech- 
niques to a set of Cepheid-bearing open clusters that have 
good multi-color photometry. Selection of cluster sample and 
compilation of cluster photometry are presented in § 2. Proce- 
dures on the isochrone calibration and MS-fitting method are 
described in § 3. New cluster distances and extinctions are 
presented in § 4. In § 5 we construct multi-wavelength Galac- 
tic P-L relations and discuss various systematic errors. In § 6 
we estimate distances and reddenings for the maser galaxy 
NGC 4258, the LMC, and M33. We summarize our results in 
§ 7 and discuss the importance of extinction corrections in the 
Cepheid distance scale. 

2. CLUSTER AND CEPHEID DATA 
2.1. Cluster Selection, Metallicity, and Age 

About 30 open clusters and associations in the Galaxy are 
known to harbor Cepheid variables (Feast & Walker 1987; 
Feast 1999), but not all of these are useful for MS fitting be- 
cause many are extremely sparse or poorly studied. We first 
excluded systems with overtone or double-mode Cepheids 
and then examined the available photoelectric and/or CCD 
photometry for each cluster. From visual inspection of color- 
magnitude diagrams (CMDs) we narrowed down the list to 
10 promising candidates for MS fitting. The populous twin 
clusters h&x Per were excluded from the current analysis be- 
cause their young age introduces significant uncertainties in 
the isochrone calibration, and the membership of UY Per in 
the clusters is also doubtful (see Walker 1987, and references 
therein). The eight remaining clusters are listed in Table 1 
with their 10 Cepheid variables. We denoted anonymous van 
den Bergh (C0634+031; van den Bergh 1957) as VDB 1 for 
brevity. Most of them are close enough that Two Micron 
All Sky Survey (2MASS) (Skrutskie et al. 2006) photometry 
reaches well down the MS. 

It is likely that all Cepheids in Table 1 are members of their 
associated clusters. For the three Cepheids in NGC 7790 there 
exist no radial velocity or proper-motion measurements, al- 
though all of them lie near the optical center of the cluster 
(Sandage 1958). V Cen is located about 25' away from the 
center of NGC 5662, but Baumgardt et al. (2000) derived a 
high membership probability from Hipparcos proper motions. 
The rest of the Cepheids are generally considered to be mem- 
bers of their clusters since they are located within the optical 
radius defined by bright early-type MS stars and red giants. In 



many cases, membership is also supported by radial velocities 
(Mermilliod et al. 1987, and references therein) and proper 
motions (Baumgardt et al. 2000, and references therein). Or- 
satti et al. (2001) suggested that TW Nor may not be a member 
of Lynga 6 from polarization measurements, but we neverthe- 
less included it here. 

Previous studies have neglected the effects of cluster metal- 
licity and age on MS fitting (e.g., Feast & Walker 1987; Hoyle 
et al. 2003). However, the luminosity of the MS at a fixed 
color (or T e ff) is sensitive to the metal abundance by A(m- 
M) /A[M/H] ~ 1 (Paper III). In addition, the cluster age 
changes the mean color of the upper MS by A(B-V) ~ 0.05 
between zero-age MS (ZAMS) and ^100 Myr isochrone, 
which is more appropriate for Cepheid-bearing clusters (see 
below). This, in turn, could affect the reddening estimates 
from the upper MS by AE(B-V) ~ 0.05 and an MS-fitting 
distance by A(m-M) <~ 0.1. 

We adopted Cepheid metal abundances from the high- 
resolution spectroscopy of Fry & Carney (1997) and assumed 
the same metallicities for the clusters. In support of this as- 
sumption, Fry & Carney derived metallicities for two dwarfs 
in M25 that are in agreement with the metallicity of U Sgr. 
In addition, Fry & Carney measured ([Fe/H]) = +0.01 ± 0.02 
for two Pleiades dwarfs, which is in good agreement with our 
adopted value for the Pleiades, [Fe/H] = +0.04 ± 0.02. 

Fry & Carney also derived a-element abundances (Si, Ca, 
and Ti) for their sample. For the stars in this paper, the aver- 
age enhancement is ([a/Fe]) = 0.14 ±0.02. We determined 
an effective metallicity [M/H] from the measured [Fe/H] and 
[a/Fe] using the Kim et al. (2002) procedure. These abun- 
dances are shown in the fourth column of Table 1 . The mean 
metal abundance of our sample is ([M/H]) = +0.03 ±0.03 
and the standard deviation is 0.08 dex, indicating an intrinsic 
spread in metallicity for our clusters. On average, the rescaled 
[M/H] are larger than [Fe/H] by 0.06 dex. We computed er- 
rors in [M/H] by propagation of errors in [Fe/H] and [a/Fe]. 

The metallicity scale for the Cepheids is probably robust to 
^0.1 dex. As part of a larger study, Yong et al. (2006) rean- 
alyzed spectra of 1 1 Cepheids in Fry & Carney and derived 
lower [Fe/H] by m 0. 1 dex and lower a-element enhancement 
by 0.05 dex on average. The mean difference in the to- 
tal metallicity is A([M/H]) = -0.15 ± 0.03, in the sense that 
the new analysis indicates lower values. Independently, An- 
drievsky et al. (2002) measured abundances of 19 Cepheids in 
common with Fry & Carney. The average difference between 
these two studies is 0.1 dex in [Fe/H], the Andrievsky et 
al. values being larger. They also showed that the alpha abun- 
dance ratio ([Si, Ca, Ti/Fe]) is nearly solar, which results in a 
small difference in the total metal abundance between Fry & 
Carney and Andrievsky et al., A ([M/H]) = -0.01 ±0.02. 

We adopted a single representative age (80 Myr) for 
all clusters, based on comparing our isochrones (§ 3.1) 
to the cluster CMDs (§ 2.2), and took its uncertainty as 
A log t (Myr) = 0.2 (f = 80^q Myr). The uncertainty represents 
both a range of ages among the clusters and the accuracy of 
the fitting. The average age in the Lynga Open Clusters Cat- 
alog (Lynga 1987) is 72 Myr for our sample clusters with a 
standard deviation of 38 Myr. Previous age estimates based on 
isochrone fitting are generally within our adopted age range. 
Mermilliod (1981) found that NGC 7790 and NGC 6067 are 
about as old as the Pleiades (78 Myr) from models with core 
overshooting. Meynet et al. (1993) estimated w 70 Myr for 
NGC 5662 and NGC 6087 when the Pleiades is 100 Myr old 
from core-overshooting models, while the age of NGC 6067 
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TABLE 1 

List of Cepheid Clusters and Cepheid Properties 



Cluster 


Other Name 11 


Cepheid 


[M/H] b 


logP 


(B) 


(V) 


ik) 


(J) 


(H) 


(K s ) 


NGC 129 


C0027+599 


DL Cas 


+0.08 ±0.03 


0.903 


10.119 


8.969 


7.655 








VDB l c 


C0634+031 


CV Mon 


+0.07 ±0.08 


0.731 


11.607 


10.304 


8.646 


7.310 


6.783 


6.475 


NGC 5662 


C1431-563 


VCen 


+0.00 ± 0.03 


0.740 


7.694 


6.820 


5.805 


5.006 


4.623 


4.462 


Lynga 6 


C1601-517 


TW Nor 


+0.12 ±0.05 


1.033 


13.672 


11.667 


9.287 


7.403 


6.686 


6.225 


NGC 6067 


CI 609-540 


V340 Nor 


-0.14 ±0.05 


1.053 


9.526 


8.370 


7.168 


6.192 


5.724 


5.515 


NGC 6087 


C1614-577 


SNor 


+0.05 ±0.03 


0.989 


7.373 


6.429 


5.422 


4.661 


4.269 


4.115 


M25 


C1828-192 


USgr 


+0.10 ±0.06 


0.829 


7.792 


6.695 


5.448 


4.511 


4.085 


3.892 


NGC 7790 


C2355+609 


CEa Cas 




0.711 


12.070 


10.920 


9.470 d 








NGC 7790 


C2355+609 


CEb Cas 




0.651 


12.220 


11.050 


9.690 d 








NGC 7790 


C2355+609 


CF Cas 


-0.11 ±0.05 


0.688 


12.335 


11.136 


9.754 









a "C" cluster designation from Lynga (1987). 

b An effective metallicity estimated from [Fe/H] and [ct/Fe] in Fry & Carney (1997) using the Kim et al. (2002) procedure. 
c Also named CV Mon cluster or anonymous VDB cluster (van den Bergh 1957). 
d Taken from Tammann et al. (2003). 



was estimated to be w 170 Myr. Cluster ages can also be 
inferred from the masses of Cepheids. There are several clas- 
sical Cepheids in binary systems, and these have dynamical 
masses of about 5-6M Q (Bohm-Vitense et al. 1998; Evans & 
Bolton 1990; Evans et al. 1997, 1998, 2006) over the period 
range similar to those of our sample [0.65 < logP(days) < 
1.05]. For our models, stars of this mass range have an MS 
lifetime of about 50-80 Myr. 

2.2. Photometry 

For the photometry of the Cepheids, we used intensity- 
mean apparent magnitudes in BVIc from Berdnikov et al. 
(2000), which are on the Cape (Cousins) system as realized 
by Landolt (1983, 1992). For CEa Cas and CEb Cas, the I c - 
band measurements were taken from Tammann et al. (2003). 
We adopted JHK photometry from Laney & Stobie (1992) on 
the Carter (1990) system and transformed it to the LCO sys- 
tem using the Persson et al. (2004) transformation equations. 
We assumed 0.025 mag zero-point error in the Cepheid pho- 
tometry and treated this error independently from the cluster 
photometry (see below). Table 1 lists the Cepheid photome- 
try. 

For the cluster photometry, we compiled photoelectric and 
CCD BVIc data on the Johnson-Cousins system from the lit- 
erature, as well as from WEBDA (Mermilliod & Paunzen 
2003). 1 The references for the photometry are shown in Ta- 
ble 2, where we used WEBDA's cross-identification of opti- 
cal sources among different references except for a few cases 
where we found missing entries and misidentifications. For 
each cluster we picked one or two references to define a local 
standard in V or B - V by weighing factors such as the num- 
ber of observations, number of stars with photometry, mag- 
nitude range in V, the photometric calibration procedure, and 
whether the photometry was generally in agreement with data 
in other studies. The remaining columns of Table 2 show 
the mean differences in V and B - V with respect to the lo- 
cal standard, and the number of stars in common. The errors 
shown are standard errors of the mean difference. For VDB 1 
we combined the photoelectric photometry by Arp (1960) and 
CCD photometry by Turner et al. (1998, their Table 2), whose 
photometry was tied to the former study. 

In most cases the differences between one study and an- 

1 See http://obswww.unige.ch/webda/webda.html. 



other showed no trends with magnitude. In M25, however, 
the differences in B-V between Johnson (1960) and Sandage 
(1960) were statistically significant (> 2.5cr): 

A(B-V)(Sandage-Johnson) = +0.009(V-10)-0.010. 

Similarly, in NGC 6087 there were significant trends with 
magnitude for Fernie (1961) and Breger (1966) with respect 
to Turner (1986): 

AV (Fernie- Turner) = -0.033(V- 10)-0.052, 

AV (Breger -Turner) = -0.014(y- 10) + 0.002. 

We applied these corrections to put them on the same scale 
as the reference photometry. For NGC 7790 we only used 
the photometry by Stetson (2000), although we computed the 
differences with respect to other studies and included these in 
Table 2. We applied zero-point corrections to individual val- 
ues if the average difference was significant at the 2.5cr level. 
After shifting to a common scale, magnitudes and colors in 
multiple studies were averaged together. The weighted rms 
differences between the studies in Table 2 are 0.026 mag in V 
and 0.024 mag in B-V. We therefore adopted 0.025 mag as 
the characteristic size of systematic errors in the photometry. 

Useful V -lc photometry for MS fitting is limited only to 
four clusters: VDB 1 (Turner et al. 1998), Lynga 6 (Walker 
1985), NGC 6067 (Piatti et al. 1998), and NGC 7790 (Romeo 
et al. 1989; Gupta et al. 2000; Stetson 2000). As with B-V we 
adopted the Stetson (2000) photometry for NGC 7790. Com- 
pared to Stetson's data, the photometry in Romeo et al. (1989) 
is bluer by 0.026 ± 0.003 mag, and the photometry in Gupta et 
al. (2000) is redder by 0.023 ± 0.001 . Assuming these values 
as a characteristic size of error, we adopted 0.025 mag for the 
photometric zero-point error in V -lc- 

We combined optical photometry with JHK S measurements 
from the All Sky Data Release of the 2MASS Point Source 
Catalog (PSC) 2 . WEBDA provides celestial coordinates for 
only a small fraction of stars in this study. For the others 
we identified each source on the images from the Digitized 
Sky Survey, 3 or we computed the celestial coordinates from 

2 See http://www.ipac.caltech.edu/2mass/. 

3 The Digitized Sky Surveys were produced at the Space Telescope Sci- 
ence Institute under US Government grant NAG W-2166. The images of 
these surveys are based on photographic data obtained using the Oschin 
Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. 
The plates were processed into the present compressed digital form with the 
permission of these institutions. 
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TABLE 2 

Photometry References and Zero-Point Differences 



Reference AV 






A(B-V) 


Ncomp 


NGC 129 


Arp et al. (1959) 


+0.010 ±0.006 


26 


+0.004 ±0.011 


26 


Hoag et al. (1961) 


+0.000 ±0.011 


18 


+0.059 ±0.008 


18 


Turner et al. (1992) 


+0.015 ±0.006 


20 


+0.048 ±0.008 


22 


Phelps & Janes (1994) 


standard 




standard 




VDB 1 


Arp (1960) 


standard 




standard 




Turner et al. (1998) 


tied to Arp 




tied to Arp 




NGC 5662 


Moffat &Vogt (1973) 


+0.009 ± 0.005 


13 


+0.002 ± 0.004 


28 


Haug (1978) 


standard 




+0.027 ± 0.004 


20 


Claria et al. (1991) 


+0.051 ±0.004 


19 


standard 




Sagar& Cannon (1997) 


-0.021 ±0.015 


10 


+0.004 ±0.006 


14 


Lynga 6 


Madore (1975) 


+0.009 ±0.015 


8 


-0.002 ±0.010 


9 


Moffat &Vogt (1975) 


-0.097 ±0.018 


9 


-0.030 ±0.025 


10 


van den Bergh & Harris (1976) 


-0.010 ±0.008 


12 


+0.024 ±0.014 


12 


Walker (1985) 


standard 




standard 




NGC 6067 


Thackeray et al. (1962) 


-0.019 ±0.013 


15 


+0.012 ±0.008 


14 


Walker & Coulson (1985) 


standard 




standard 




Piattietal. (1998) 


-0.021 ±0.003 a 


123 






NGC 6087 


Fernie (1961) 


-0.057 ±0.020 


9 


+0.041 ±0.006 


10 


Breger(1966) 


+0.005 ±0.008 


18 


-0.019 ±0.011 


17 


Turner (1986) 


standard 




standard 




Sagar& Cannon (1997) 


+0.080 ± 0.070 


2 


+0.025 ± 0.035 


2 


M25 


Johnson (1960) 


-0.018 ±0.003 


74 


standard 




Sandage (1960) 


+0.018 ±0.005 


51 


+0.005 ± 0.004 


68 


Wampler et al. (1961) 


standard 




+0.021 ±0.004 


51 


NGC 7790 


Sandage (1958) 


-0.010 ±0.013 


26 


-0.018 ±0.008 


26 


Romeo et al. (1989) 


-0.001 ±0.003 


186 


+0.013 ±0.003 


184 


Phelps & Janes (1994) 


-0.030 ± 0.001 


360 


-0.005 ± 0.002 


337 


Lee & Lee (1999) 


+0.027 ± 0.003 


164 


-0.032 ± 0.003 


157 


Gupta et al. (2000) 


+0.008 ±0.001 


192 


+0.009 ±0.002 


187 


Stetson (2000) 


standard 




standard 





NOTE. — The mean color and magnitude differences are computed after 3<7 rejection. Uncertainties 
are the standard error of the mean. The differences are in the sense of individual values minus those 
adopted as the local standard. The photometry in Turner et al. (1998) are directly tied to the photometry 
of Arp (1960). 

"Comparison for stars with V < 15. 



the plate position information in WEB DA. The rms difference 
between the retrieved and the 2MASS coordinates was typi- 
cally 0.5". For NGC 6067 we matched 2MASS sources with 
optical photometry only for those stars with good positional 
accuracy. The validation of the 2MASS source matches was 
confirmed from the resulting tight optical and near-infrared 
color-color relations. Based on PSC flag parameters, we ig- 
nored the infrared data if sources were undetected, blended, 
or contaminated. Calibration errors in JHK S were taken as 
the uncertainty specified in the explanatory supplement to the 
2MASS All Sky Data Release: 0.01 1 mag in J, 0.007 mag in 
H, and 0.007 mag in K s . 4 

Most studies do not report individual errors in optical mag- 
nitude or color, so we computed the median of V and B - V 
differences from various studies (after shifting to a common 
scale using the > 2.5a criterion) and assigned this value to 
represent the random photometric errors for all stars in each 
cluster. For V— Iq we assigned 0.02 mag as the random er- 
ror. Errors for individual stars in Stetson (2000) were adopted 

4 See http://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html. 



without change. The V — K s errors were computed as the 
quadrature sum of V errors and the catalog's "total" photo- 
metric uncertainties in K s . 

The cluster CMDs in various color combinations are shown 
in Figures 1 through 4. Stars rejected as being far from the MS 
are shown as open circles, while those used in the MS fitting 
are shown as bull's-eyes. The curved line in each panel is the 
best-fitting isochrone on each CMD. Details of the MS fitting 
and outlier rejection are given in the next section. Through- 
out this paper we denote CMDs by their color and luminosity 
indices as (B-V,V), (V-I C ,V), (V-K S ,V), (J-K S ,K S ), and 
(H-K S ,K S ). 

3. MAIN-SEQUENCE FITTING 

3.1. Extension of Isochrone Calibration 

In our first two papers of this series (Paper I and Paper II), 
we assessed the accuracy of distances from MS fitting and 
examined systematic errors in the transformation of theoret- 
ical to observational quantities. In Paper I we demonstrated 
that stellar models from the Yale Rotating Evolutionary Code 
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FIG. 1 . — Cluster CMDs in B— V. The bull's-eyes are stars remaining after the photometric filtering, while open circles are stars rejected. The solid lines are 
empirically calibrated isochrones at the best-fitting distance (Tables 6), reddening (Tables 7), and R v q (Tables 9). The boxes on the upper MS and lower MS are 
regions where the color excess and distance were determined, respectively. 




(Sills et al. 2000) are in agreement with the masses and lu- 
minosities of the well-studied Hyades eclipsing binary vB 22 
(Torres & Ribas 2002). These models also satisfy stringent 
tests from helioseismology and predict solar neutrino fluxes 
in line with observations (Basu et al. 2000; Bahcall et al. 
2001; Bahcall & Pinsonneault 2004). In Paper II we showed 
that the models provide a good match to the spectroscopi- 
cally determined temperatures (Paulson et al. 2003) of in- 
dividual Hyades members with good parallaxes (de Bruijne 



et al. 2001). However, we found that any of the widely-used 
color- T e ff relations (e.g., Alonso et al. 1995, 1996; Lejeune et 
al. 1997, 1998, hereafter LCB) fail to reproduce the observed 
shapes of the MS in the Hyades, having differences in broad- 
band colors as large as 0.1 mag. The existence of these sys- 
tematic errors in the colors, in the presence of agreement be- 
tween the spectroscopic and theoretical luminosity-r e ff scales, 
led us to argue that the problems lie with the adopted color- 
T e g relations instead of errors in the theoretical T e ff scale. We 
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FIG. 3. — Same as Fig. 1, but CMDs in V-K s . The dotted line represents the 2MASS completeness limit. 



then adopted the LCB color- r eff relations and computed em- 
pirical corrections to the LCB relations to match photometry 
in the Hyades. In Paper III we showed that isochrones with 
the empirical color- T e ff corrections accurately match the MS 
shapes of other nearby clusters in three color indices (B - V, 
V-Ic, and V-K s ), yielding distance estimates with errors as 
small as 0.04 mag in distance modulus (or 2% in distance). 

However, the Hyades-based calibration is of limited use for 
this study. The calibration is reliable only at My > 3 because 
the number of stars in the upper MS of the Hyades is small in 
this relatively old open cluster (550 Myr from models exclud- 
ing convective overshoot; Perryman et al. 1998). On the other 
hand, photometry of our sample clusters is quite sparse or ab- 
sent in this lower MS part. It is therefore desirable to define 
an extension of the empirical color- T e ff corrections for hot- 
ter and brighter stars. A good choice for a calibrating cluster 
is the Pleiades for several reasons. First of all, its geometric 
distance is now known to a precision comparable to that of 
the Hyades, (m-M) = 5.63 ± 0.02 (see Paper III and refer- 
ences therein for a detailed discussion of the cluster parame- 
ters). Second, it has low reddening, E(B-V) = 0.032 ±0.003. 
Third, a precise metal abundance from high-resolution spec- 
tra is available, [Fe/H] = +0.04 ±0.02. Furthermore, Cepheid 
clusters tend to be of ages comparable to that of the Pleiades, 
making it a good template for distance and reddening esti- 
mates. 

The sources of our Pleiades photometry are described in Pa- 
per III. We began by removing a few known binaries and then 
rejected stars that are more than 0.1 mag in color away from a 
15 point median of the sample sorted in V; this step was inde- 



pendent of the theoretical isochrones. Table 3 lists photome- 
try of the hot Pleiades stars that remained after this selection. 
The CMDs of these stars are shown in Figure 5. Small plus 
signs denote stars rejected as far from the MS, while open cir- 
cles show stars retained for use in the following calibration 
process. 

In the course of this work we became aware of some prob- 
lematic features in the LCB color tables. They computed 
synthetic colors from both the "original" and the "corrected" 
model flux distributions (hereafter original and corrected LCB 
tables, respectively). In Figure 6 we compare the corrected 
LCB table with the original LCB table for MS stars (logg = 
4.5). The LCB empirical color- r e ff relations (corrected LCB 
table) are defined only for r e ff < 1 1500 K, which results in an 
artificial jump in CMDs, particularly in the V-band bolomet- 
ric correction. In addition, there are small scale structures in 
the color corrections (especially in V—Ic) that cause interpo- 
lation noise in our isochrones at T e ff > 8000 K. We therefore 
used the original LCB table above 8000 K, the corrected LCB 
table below 6500 K, and a linear ramp between the two tables 
in order to produce smoother base isochrones. Isochrones at 
an age of 100 Myr generated with these "merged" color- T e ff 
relations are shown as dot-dashed lines in Figure 5. However, 
these still fail to match the observed MS of the Pleiades and 
therefore require further corrections. 

We derived empirical corrections to the isochrone in B - V , 
V—Ic, V — K s , and J-K s by forcing a match to the Pleiades 
photometry as done for the Hyades (Paper II). Figure 7 
plots the difference in color between individual stars and an 
isochrone generated with the merged LCB color- r e ff relations. 
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The r eff was calculated from M v . Filled circles are for the 
Pleiades stars, and the open circles are for the Hyades stars 
used in the calibration of the lower MS. The Pleiades stars 
have smaller error bars than the Hyades points because we 
assumed a negligible distance spread for the Pleiades, while 
the Hyades ones reflect parallax errors for individual stars (de 
Bruijne et al. 2001). The solid lines in Figure 7 represent 
our empirical corrections, which were computed at a con- 
stant My in a seven-point moving window. Given the modest 
sample size and intrinsic scatter in the upper MS, we did not 
attempt to resolve out smaller scale structures. We defined 
the empirical color- T eB corrections as those for the Pleiades at 
r eff > 6838 K and for the Hyades at T eB < 6500 K. A spline 
function was used to bridge over the narrow r eff gap between 
these two corrections. We also used H-K s in MS fitting but 
without color corrections because photometric errors were too 
large in this color to provide a reliable calibration. 

Applying these corrections to the colors at each T e ff de- 
fines the empirically calibrated isochrone for the Pleiades, 
which is tabulated along with the color corrections in Table 4. 
Isochrones incorporating the corrections for the Pleiades are 
shown as solid lines in Figure 5. Systematic errors in the color 
corrections are ~ 0.01 mag from the errors in the adopted dis- 
tance, metallicity, and reddening values for the Pleiades. We 
assumed that our empirical corrections in the upper MS are 
independent of metallicity and age and applied them to all 
isochrones generated in Paper III. The isochrones constructed 
in this way are available online at 

http : // www . astronomy .ohio-state.edu/iso/. 
There is one important functional difference between the 



calibration for luminous stars and that for the lower MS stars. 
The MS for FGK-type stars is intrinsically very narrow. As 
a result, for the low-mass stars we could plausibly claim that 
the empirical cluster sequence could be used to redefine the 
color- T e ff relationship. On the other hand, the dispersion about 
the mean trend for upper MS stars is significantly larger. For 
the Pleiades, we found a dispersion of 0.02 mag in B - V for 
0.0 < B — V < 0.40 (after excluding outliers), which is sig- 
nificantly above the photometric precision. A comparable in- 
crease was seen in V— Ic and V — K s (0.02 and 0.04 mag, re- 
spectively). Although this could be due to differential redden- 
ing across the cluster field (Breger 1986), the persistence of 
rapid rotation in early-type stars can impact both their evo- 
lution and the mapping of T e g onto colors (e.g., Collins & 
Sonneborn 1977) as a consequence of the Kraft (1965) break 
in rotational properties. Furthermore, these effects depend on 
both the rotation rate and the inclination to the line of sight 
and can be as large as 0. 1 mag in B-V . Nonetheless, the mean 
cluster locus can still be used to define a template for distance 
studies but should not be interpreted as a direct change in 
the color- T e ff relationship. In other words, our empirical cor- 
rections remove systematic trends from the rotation-induced 
color anomalies, as well as from any intrinsic problems in the 
adopted color- T e g relationship. 

3.2. Main-Sequence Fitting 

We determined the distance, reddening, and the ratio of to- 
tal to selective extinction in V [Ry = Av/E(B-V)] for each 
cluster via an iterative approach, which fits an isochrone to 
individual CMDs of various color indices. All of these pa- 
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TABLE 3 

Merged Photometry in the Pleiades 
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1.002 
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-0.027 ±0.021 
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1.377 
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0.010 
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"The star designation is that of Hertzsprung (1947), Hz II; and Trumpler (1921), TS. 



TABLE 4 

Empirically Calibrated Pleiades Isochrone 



Teff 

(K) 


% 


(B-V)o 
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(J-K s ) 


A(B-V) 
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A(V-K S ) 


A(J-K S ) 
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0.028 
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0.005 


7853 


2.487 


0.221 


0.269 


0.527 


0.084 


-0.006 


0.030 


0.037 


0.002 


7606 


2.652 
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0.925 
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0.015 
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FIG. 5. — Empirically calibrated isochrones for the Pleiades (solid lines). The dot-dashed lines show the isochrones without the color corrections. The open 
circles are stars retained for use in the calibration process, and the small plus signs show stars rejected (see text). 
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FIG. 6. — Differences in the bolometric correction and colors between 
"corrected" and "original" LCB tables for MS stars ([M/H] = 0, logg = 4.5). 
Differences are in the sense of the former minus the latter values. 

rameters are essential to obtain self-consistent solutions for 
the absolute magnitude of a Cepheid (§ 5). Our MS-fitting 
process also includes rejection of stars that are too far away 
from an MS ridge-line, defined by the locus of points with 
the highest density on each CMD. Individual color excesses 
in each cluster were used to derive R v according to the extinc- 
tion law described by Cardelli et al. (1989, hereafter CCM89). 
The rest of this section describes the MS-fitting process in de- 
tail. 

Our sample clusters do not have detailed radial velocity or 
proper-motion membership studies. Before isochrone fitting 
we therefore applied the photometric filtering procedure de- 
scribed in Paper III to identify and reject stars that are likely 
foreground/background objects or cluster binaries. The pro- 
cedure iteratively identifies a cluster MS ridge-line indepen- 
dently of the isochrones and rejects stars if they are too far 
away from the ridge-line. This was done by computing x 2 in 
each CMD as 



X 



1=1 



N 

£ 

1=1 



(AX;) 2 



4, ! + (7, 1 c r v,«) 2 + CTo 



2 ' 



(1) 



where AX, is the color difference between the i' h data point 
and MS at the same V; ox.i and av.i are photometric errors in 
color and V, respectively. The error in V contributes to the 
error in AXj by the inverse slope of the MS, 7^. We added 
do in quadrature to the propagated photometric errors in the 
denominator to take into account the presence of differential 
reddening, cluster binaries, non-cluster members, and other 
effects that would increase the MS width above the photomet- 
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FIG. 7. — Empirical color- T c ff corrections (solid line). The filled circles 
show the differences in colors between the Pleiades data and the uncorrected 
isochrone at a constant My . The sense of the difference is the former minus 
the latter values. The open circles are those for the Hyades stars. 

ric precision. We adjusted the value of cto so that the total x 2 
is equal to the total number of stars N used in the filtering. 

Initially we rejected all data points as outliers if xj (the in- 
dividual contribution to x 2 ) was greater than 9 (correspond- 
ing to a 3er outlier). We repeated adjusting ctq and reject- 
ing outliers with the reduced set of data points until there re- 
mained no point with x 2 > 9. We combined the results from 
all CMDs, and rejected stars if they were tagged as an out- 
lier at least in one of the CMDs. Missing data in any colors 
was not a criterion for rejection. For this study we included 
(J-K S ,K S ) and (H-K S ,K S ) in the filtering procedure. We 
reduced the rejection threshold at each iteration until it was 
limited to 2.0-2.4cr, but our result is insensitive to the final 
rejection threshold within the fitting errors. 

Our filtering results are shown in Figures 1-4. The bull's- 
eyes represent stars retained and used in the following MS fits, 
and open circles are the ones rejected by the filtering. At the 
end of the iteration, <jq values in (B — V,V) and (V — K S ,V) 
were (0.055 mag, 0.101 mag) in NGC 129, (0.050 mag, 
0.102 mag) in VDB 1, (0.023 mag, 0.064 mag) in NGC 5662, 
(0.048 mag, 0.102 mag) in Lynga 6, (0.031 mag, 0.078 mag) 
in NGC 6067, (0.033 mag, 0.058 mag) in NGC 6087, 
(0.043 mag, 0. 1 1 1 mag) in M25, and (0.01 1 mag, 0.026 mag) 
in NGC 7790. These excess dispersions are most likely due to 
differential reddening of heavily reddened, young open clus- 
ters. Note that the ratio of cr values from (B-V,V) and 
(V—K S ,V) is also approximately equal to the color-excess ra- 
tio E(V-K S )/E(B-V) w 2.8 (see below). We discuss in detail 
the effects of differential reddening on P-L relations in § 5. In 
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(J—K S ,K S ) and (H — K S ,K S ), photometric errors were usually 
large enough that <tq values were automatically set to zero. 

The fitting process took advantage of the changing slope 
of the MS with absolute magnitude. Near the MS turnoff 
(My < 1 .5), the MS is nearly vertical, allowing a precise de- 
termination of the reddening. Further down (1 .5 < My < 3.5), 
the slope of the MS is still relatively steep [e.g., AMy / A(B - 
V) ~ 5], but it nevertheless provides sufficient leverage for 
determining the distance. We divided each CMD into two 
zones, shown as boxes in Figures 1-4, which define the re- 
gions used for determining the color excess and distance. The 
color excess was determined for stars with -1 .0 < My < 2.0 in 
(fi-V,V), (V-/ c ,V),and (V-K S ,V), and at -0.8 <M K < 1.7 
in (J-K S ,K S ) and (H-K S ,K S ). Distance determinations came 
from stars with -0.1 < (B-V)q < 0.5 and the corresponding 
color ranges in other colors with the same My. For Lynga 6 
and NGC 7790, (V-K S ,V) and (J-K S ,K S ) were not used in 
the distance estimation because the lower parts of the MS are 
below the 2MASS completeness limit (K s w 14.3). 

Since the color excess and R v depend on the intrinsic 
color of stars, we adopted the color-dependent relations from 
Bessell et al. (1998, hereafter BCP98) and M. S. Bessell 
(2007, private communication). Their formulae were based 
on the theoretical stellar spectral energy distributions with the 
extinction law from Mathis (1990). Specifically, their formu- 
lation yields 

R v = R vfi + 0.22(B- V) 0) (2) 



E(B -V) = E(B - V) { 1 - 0.083(5 - V) }, 



(3) 



where (B - V)o is an intrinsic color of each star. The /?y,o and 
E(B - V)o are the values for zero-color stars, where Ry o is 
3.26 in their model. For other colors we used the following 
color-excess ratios: 



E(V- 


-Ic) 


E(B- 


-V) 


E(V- 


-K s ) 


E(B- 


-V) 


E(J- 


-K s ) 


E(B 


-V) 


E(H- 


-K s ) 


E(B- 


-V) 



1.30 + 0.06(V- 


-/c)o 


(4) 


2.88 + 0.07(V- 


-K s ) 


(5) 


0.56 + 0.06(7- 


K s )o 


(6) 


0.20+0.15(ff- 


~K s )q. 


(7) 



The first terms represent the color-excess ratios from CCM89 
at Ryfi = 3.26, and the second terms are color-dependent re- 
lations from BCP98. Differences between the original coeffi- 
cients in BCP98 and those for the 2MASS JHK S system using 
transformation equations in Carpenter (200 1) 5 were found to 
be negligible. Table 5 lists color-excess ratios for zero-color 
stars from CCM89 and BCP98, as well as observationally 
derived values from Schultz & Wiemer (1975), Dean et al. 
(1978), Sneden et al. (1978), and Rieke & Lebofsky (1985). 
Color-excess ratios determined from our cluster sample are 
also shown in the fourth column, which is discussed in the 
next section. 

The fitting process begins with a guess at the average dis- 
tance. At this distance, individual color excesses for zero- 
color stars [£(fi-V)° bs , E(V - / c )g bs , etc.] were determined 
using stars on the upper MS. These were used to determine 

5 Updated color transformations for 2MASS All-sky Data Release can be 
found at 

http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_4b.html. 



the best-fit R vfi and E(B - 
law by minimizing 



V)o from the CCM89 extinction 



X 



£(A,-y)S bs -£(A,-y)^ CM89 



0~E(\-V) 



(8) 



where A, represents effective wavelengths in BIq (Bessell 
et al. 1998) and isophotal wavelengths in JHK S (Cohen et al. 
2003). The index in the summation runs over B, Iq, J, H, and 
K s . The (Je(x-v) is the fitting error in the upper MS, which was 
determined from the scatter of the points around the isochrone 
and the individual errors in color. The CCM89 color excess, 
E(X- V)g CM89 , can be expressed as 



E(X- 



T/ xCCM89 



■ Ry. Q E(B-V) 



Ax 
Ay' 



(9) 



where A \ /Ay is a function of Ryp from polynomial relations 
in CCM89. Note that R vfi and E(B-V)q in equation (8) are 
two parameters to be determined and that E(B-V)o in equa- 
tion (9) is the normalization factor, which is different from the 
input color excess in (B - V, V). 

The fitting procedure then found weighted median distances 
on each CMD by using the stars in the lower part of the MS. 
The weighted average distance from all CMDs was computed, 
and the process was repeated until convergence, which oc- 
curred when the difference in the average distance from the 
previous iteration became smaller than its propagated error. 

An example of the fit is shown in Figure 8 for NGC 6067. 
The solid line in the top panel is the best-fitting CCM89 ex- 
tinction curve to the data (open circles). The bottom panel 
shows likelihood contours in R v $ and E(B-V)o drawn at A% 2 
= 2.30, 6.17, and 11.8 from the minimum \ 2 = 4.5 (68.3%, 
95.4%, and 99.73% confidence levels for 2 degrees of free- 
dom). 

4. CLUSTER PROPERTIES 

In this section we present our results on the distance, red- 
dening, and Ry of the sample clusters and evaluate various 
systematic errors in these parameters. In particular, we assess 
the accuracy of MS-fitting based on internal consistency of 
distances from multi-color CMDs. We show that our distance 
estimates are generally shorter than previous ZAMS-fitting re- 
sults and that this is mainly due to our lower reddening esti- 
mates. 

4.1. Distance, Reddening, and Ry 

Our MS-fitting results on the distance and color excess from 
individual CMDs are shown in Tables 6 and 7, respectively, 
along with a number of stars used in the fit. The last two 
columns in Table 6 show the weighted average distance mod- 
ulus from all available CMDs with its propagated error and 
the standard deviation (a) of individual distances. The last 
column in Table 7 lists a standard deviation of E(B - V) from 
all CMDs, after transforming individual color-excess values 
in V-I c , V-K s , and J-K s to E(B-V) using the CCM89 
reddening law (Table 5). We excluded E(H-K S ) from this 
computation since we found a large systematic deviation of 
our color-excess estimates from the CCM89 law as discussed 
below. 

In Figures 1-4 the best-fitting isochrones are overlaid on 
the cluster CMDs. While there is some variation in the qual- 
ity of the data, we could divide the sample clusters into three 
groups based on how tight the MS is on optical CMDs. The 
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TABLE 5 

Color-Excess Ratios for Zero-Color Stars 



Color-Excess Ratios 


CCM89" 


BCP98" 


This Paper b 


Others c 


ref 


E(V-I C )/E(B-V) 


1.30 


1.32 


1.26 ±0.06 


1.25 ±0.01 


1 


E(V-K S )/E(B-V) 


2.88 


2.91 


2.82 ±0.06 


2.78 ±0.02 


2,3,4 


E(J-K S )/E(B-V) 


0.56 


0.58 


0.53 ±0.04 


0.52 ±0.03 


4 


E(H-K S )/E(B-V) 


0.20 


0.22 


0.16 ±0.02 


0.20 ±0.04 


4 



REFERENCES. — (1) Dean et al. 1978; (2) Schultz & Wiemer 1975; (3) Sneden et al. 1978; (4) Rieke 
&Lebofsky 1985. 



NOTE. — Color excesses involving 2MASS JHK, bandpasses are corrected for the difference in filter 
effective wavelengths. 

"Estimated at R v0 = 3.26. 
'Weighted mean and rms deviation. 

c Revised for zero-color stars using color-dependent reddening relations in BCP98. 



TABLE 6 

MS-Fitting Distance 



Cluster 








(mi- 


-M) 








(m-M) , ave a 


0"(m-M) b 


(B-V.V) 


Nu 


(V-lc.V) 




(V-K S ,V) 


Nm 




Nm 


NGC 129 . 


11.030 ±0.041 


79 






11.078 ±0.044 


58 


11.007 ±0.037 


53 


11.035 ±0.023 


0.027 


VDB 1 . . . 


11.067 ±0.097 


20 


10.744 ±0.033 


13 


10.576 ±0.062 


13 


10.706 ±0.082 


18 


10.734 ±0.026 


0.209 


NGC 5662 


9.351 ±0.025 


52 






9.290 ±0.037 


46 


9.201 ±0.043 


37 


9.307 ±0.019 


0.075 


Lynga 6 . . . 


11.604±0.114 


21 


11.415 ±0.117 


21 










11.512±0.082 


0.134 


NGC 6067 


10.991 ±0.019 


141 


11. 125 ±0.026 


120 


11.004 ±0.033 


43 


10.996 ±0.059 


35 


11.029 ±0.013 


0.064 


NGC 6087 


9.708 ±0.048 


22 






9.601 ±0.045 


20 


9.649 ±0.049 


19 


9.650 ±0.027 


0.054 


M25 


8.961 ±0.030 


45 






8.899 ±0.054 


36 


8.885 ±0.052 


35 


8.934 ±0.023 


0.040 


NGC 7790 


12.474 ±0.014 


63 


12.427 ±0.019 


57 










12.457 ±0.011 


0.033 



NOTE. — Formal fitting errors are shown for individual distance moduli with a number of points used in the fit. 

"Weighted average distance and its propagated error, not including systematic uncertainties. 
'Standard deviation of individual distance moduli. 



TABLE 7 
MS-Fitting Color Excess 



Cluster 


E(B-V) 


N m 


E(V-I c )o 


Nm 


E(V-K S ) 


Nm 


E(J-K S ) 


Nm 


E(H-K,) 


Nm 




NGC 129 . 


0.500 ±0.008 


53 






1.431 ±0.020 


45 


0.280 ±0.006 


43 


0.069 ±0.006 


41 


0.002 


VDB 1 . . . 


0.573 ±0.014 


12 


0.611 ±0.023 


9 


1.458 ±0.031 


11 


0.358 ±0.013 


10 


0.066 ±0.012 


10 


0.075 


NGC 5662 


0.277 ± 0.003 


47 






0.795 ±0.012 


41 


0.172 ±0.006 


38 


0.050 ±0.005 


37 


0.018 


Lynga 6 . . . 


1.326 ±0.012 


23 


1.698±0.017 


22 


3.659 ±0.056 


9 


0.635 ±0.020 


8 


0.200 ±0.016 


9 


0.087 


NGC 6067 


0.340 ± 0.003 


107 


0.428 ±0.007 


67 


0.935 ±0.014 


31 


0.166 ±0.008 


26 


0.025 ±0.008 


22 


0.019 


NGC 6087 


0.132 ±0.007 


30 






0.410±0.010 


25 


0.084 ±0.009 


24 


0.019 ±0.007 


24 


0.009 


M25 


0.424 ±0.009 


60 






1.160 ±0.022 


52 


0.233 ±0.007 


47 


0.077 ±0.005 


47 


0.011 


NGC 7790 


0.509 ±0.002 


35 


0.634 ±0.005 


33 


1.457 ±0.021 


27 


0.248 ±0.010 


28 


0.110±0.011 


28 


0.030 



NOTE. — Color excess values are those for zero-color stars (see text). Formal fitting errors are shown for individual color excesses with a number of points used in 
the fit. 

"Standard deviation of E(B - V) estimates from all CMDs except (H - K S ,K S ), after transforming color-excess values to E(B - V) using the CCM89 color-excess ratios. 



first group contains NGC 6067 and NGC 7790, which have 
well-populated CMDs and small fitting errors on each CMD. 
They also have a well-defined peak on the distribution of 
an individual star's distance modulus as shown in Figure 9 
for NGC 7790. The second group is NGC 129, NGC 5662, 
Lynga 6, NGC 6087, and M25. These clusters have less well- 
defined peaks, either because of strong differential reddening 
(§ 5) or because of a bias in the sampling from photoelectric 
observations at faint magnitudes. The last group is VDB 1, 
which has a sparse MS and correspondingly large distance 
uncertainties. 



We can also divide the clusters into groups by examining 
the difference in derived cluster properties from different col- 
ors. Clusters where the results from multi-color CMDs have a 
high internal consistency are strong cases, while clusters with 
substantial differences indicate problems. For the above first 
and second groups of clusters except Lynga 6, the standard de- 
viations of individual distances [<7( m _A0 in Table 6] are 0.03- 
0.08 mag, which are generally consistent with the fitting error 
on individual CMDs. On the other hand, both Lynga 6 and 
VDB 1 exhibit a larger dispersion in distance [cr ( m -M) = 0.21 
and 0.13 mag, respectively], and they also have a larger stan- 
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1.20 



0.60 - 



-0.60 - 
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NGC 6067 
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i/Mm" 1 ) 




v,o 



FIG. 8.— Example of determining R va and E(B-V) for NGC 6067. 
Top: Solid line represents the best-fitting CCM89 extinction curve to color 
excesses (open circles). Bottom: Likelihood contours in R v and E(B-V)q 
shown at A\ 2 = 2.30, 6.17, and 11.8 (68.3%, 95.4%, and 99.73% confidence 
levels for 2 degrees of freedom). 

dard deviation of E(B-V) from different colors (a in Table 7). 

However, the dispersion in Lynga 6 is consistent with the 
statistical scatter expected from the small sample size. By 
contrast, VDB 1 has a statistically unreliable result that de- 
mands its exclusion from our data set. For this cluster, the 
best distance estimates inferred from different colors disagree 
by much more than their individual error estimates would im- 
ply. One possible explanation for the internal inconsistency 
is that the cluster MS is contaminated by less reddened field 
stars as shown by spectral-type study (Preston 1964). Because 
these stars happen to lie on the cluster MS, our photometric 
filtering could have misidentified cluster members. We there- 
fore exclude VDB 1 from the following analysis and leave it 
for future studies. 

Table 8 shows the error budget for NGC 6067 as an exam- 
ple. The first column displays individual sources of error, and 
the second column shows the size of the errors adopted for 
each quantity. The third through fifth columns list error con- 
tributions to Ryfl, E(B-V)o, and (m-M)o, respectively. The 
size of systematic errors in cluster metallicity, age, and pho- 
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FIG. 9. — Distribution of individual distance moduli for NGC 7790 with 
respect to the isochrones shown in Fig. 1 and 2. The open histogram is the 
distribution of all stars in the isochrone fitting range where the distance was 
determined. The hatched histogram shows only for those that passed the pho- 
tometric filtering. The vertical line indicates the weighted median distance 
modulus. 



TABLE 8 
Error Budget for NGC 6067 



Source of Error AQuantity 

[M/H] ±0.05 

log? (Myr) ±0.2 

Helium (AY) .. . ±0.010 

A(B-V) ±0.025 

A(V-I) C ±0.025 

AV ±0.025 

AJ ±0.011 

All ±0.007 

A A.'. ±0.007 

Fitting 

Total 



AR V 



AE(B-V) A(m-M) 



±0.005 


=F0.001 


±0.056 


±0.009 


=F0.013 


=F0.020 






T0.030 


±0.264 


=F0.022 


=1=0.007 


=F0.030 


=p0.002 


±0.027 


=F0.085 


±0.000 


=F0.003 


±0.009 


±0.000 


=fo.ooi 


±0.009 


±0.000 


T0.002 


±0.011 


±0.000 


=F0.005 


±0.061 


±0.004 


±0.032 


±0.286 


±0.026 


±0.079 



NOTE. — A(B-V), A(V-I)c, AV, AJ, AH, and AK S represent zero-point 
errors in the cluster photometry. 



TABLE 9 

Summary of Cluster Parameters 



Cluster 


Rv.Q 


E(B-V) 


(m-M)o 


NGC 129 


3.31 ±0.22 


0.49 ±0.03 


11.04 ±0.05 


NGC 5662 


3.30±0.38 


0.27 ±0.03 


9.31 ±0.06 


Lynga 6 


3.16±0.10 


1.32 ±0.03 


11.51 ±0.13 


NGC 6067 


3.18±0.29 


0.34 ±0.03 


11.03 ±0.08 


NGC 6087 


3.59 ±0.85 


0.13 ±0.03 


9.65 ±0.06 


M25 


3.16±0.21 


0.42 ±0.03 


8.93 ±0.08 


NGC 7790 


3.18±0.26 


0.51 ±0.02 


12.46 ±0.11 
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FIG. 10. — Mean color-excess ratios from the cluster sample. The error 
bars are from the isochrone fitting only. The solid line is a fit to the data, 
constrained to pass through the origin. The dashed line is the color-excess 
ratio from the CCM89 reddening law. 



tometry was discussed in § 2. The quantity AY represents the 
uncertainty in the helium abundance at fixed metallicity and 
comes from the consideration of the initial solar abundance 
and those of the Hyades and the Pleiades (Paper I and Pa- 
per III). The fitting errors represent the internal precision of 
the fit and were taken as the larger value of the propagated 
error or <J{ m -M) /VN in Table 6. We estimated the fitting errors 
in Rv.o and E(B-V)o from the size of the la contours on the 
A\ 2 distribution (e.g., Fig. 8). The total systematic error is 
the quadrature sum of these errors. The R V fi, E(B-V)o, and 
(m-M)o of our cluster sample are listed in Table 9 with total 
systematic errors in these quantities. 

Independently of the CCM89 reddening law, we also de- 
rived color-excess ratios based on our cluster sample. Fig- 
ure 10 shows a correlation between E(B-V) and color-excess 
values in other color indices. Individual clusters are repre- 
sented as a point, and fitting errors are shown only. The solid 
line is a linear fit to the data, constrained to pass through the 
origin. Table 5 summarizes our best-fit slopes of these lines or 
color-excess ratios. Our values are generally in good agree- 
ment with the CCM89 law and other observational estimates. 
However, our method yields a lower E(H - K s ) / E(B -V) than 
these studies, which may reflect a zero-point offset in H-K s 
of our isochrones. Note that we did not derive color correc- 
tions for H-K s . 

4.2. Comparison with Previous Studies 

The most up-to-date compilation of Cepheid parameters in 
open clusters and associations can be found in Tammann et 
al. (2003). Their distances are from the earlier compilation by 
Feast (1999). These ZAMS -fitting distances are the same as 
in Feast & Walker (1987), but include more recent distance 
estimates for NGC 129 (Turner et al. 1992) and NGC 7790 
(Romeo et al. 1989). All of our Cepheids are listed in these 
tables, and the comparisons with distances in Feast & Walker 
and Tammann et al. are shown in the upper two panels of 
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FIG. 11. — Comparison of cluster distances with those in Feast & Walker 
(1987), Tammann et al. (2003), Hoyle et al. (2003), and Gieren et al. (2005). 
The differences are in the sense of their distances minus those in this paper. 
The error bars are quadrature sum from the two studies in the lower two 
panels. Individual distances in the top two panels were revised assuming our 
Pleiades distance scale. 

Figure 11 and listed in Table 10. Both of these studies esti- 
mated distances relative to the Pleiades but assumed different 
Pleiades distance of (m-M)o = 5.57 and 5.61 mag, respec- 
tively. We revised their distances assuming our Pleiades dis- 
tance, (m-M)o = 5.63. Since they did not explicitly present 
errors, we plotted our errors only. The average differences 
in distance are A((m-M)o> = 0.09 ±0.05 and 0.12 ±0.05 
with Feast & Walker and Tammann et al., respectively. The 
sense of the differences is that our new distance estimates are 
shorter on average. The standard deviations of the differences 
are 0.14 mag for both comparisons, while the expected size 
from our error estimates is 0.07 mag. If our error estimates are 
correct, the dispersion indicates a random error of 0.12 mag 
in distance estimates from Feast & Walker and Tammann et 
al.. 

In addition to these studies, Hoyle et al. (2003) derived 
ZAMS -fitting distances to 1 1 open clusters with new photom- 
etry in UBVK. Our comparison with this study is shown in 
Figure 11 and tabulated in Table 10. The error bars repre- 
sent quadrature sums of our and their reported error estimates. 
Their distances are marginally consistent with our estimates, 
and the average difference is A((m-M)o) = 0.04±0.07, our 
distances being shorter on average. However, we did not at- 
tempt to revise their distances assuming our Pleiades distance 
scale. They used a ZAMS from Allen (1973) for the opti- 
cal data, but they also claimed that the same results were 
obtained from the ZAMS in Turner (1979) and Mermilliod 
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TABLE 10 

Comparison of Cluster Distance with Published Results 



Cluster 


This Study 


FW87" 


TSR03 a 


HST03 


G05 b 


NGC 129 . 


11.04 ±0.05 


11.30 


11.24 


10.94 ±0.14 




NGC 5662 


9.31 ±0.06 


9.19 


9.19 




9.33 ±0.06 


Lynga 6 . . . 


11.51 ±0.13 


11.48 


11.49 


11.33 ±0.18 




NGC 6067 


11.03 ±0.08 


11.19 


11.19 


11.18±0.12 


11.25±0.19 


NGC 6087 


9.65 ±0.06 


9.87 


9.87 




10.02 ±0.03 


M25 


8.93 ±0.08 


9.09 


9.09 


9.08±0.18 


8.98 ±0.02 


NGC 7790 


12.46±0.11 


12.45 


12.71 


12.58±0.14 





NOTE. — FW87, Feast & Walker (1987); TSR03, Tammann et al. (2003); HST03, Hoyle et al. 
(2003); G05, Gieren et al. (2005). 

"Distances revised assuming the same Pleiades distance, (m-M) = 5.63, as in this paper. 
b Distances to individual Cepheids from the surface brightness technique. 



(1981). Since the Pleiades distance modulus of 5.56 mag was 
adopted in Turner, the distance moduli in Hoyle et al. would 
then become longer by w 0.07 mag on average than ours when 
they are on the same Pleiades distance scale as in this paper. 

The bottom panel in Figure 1 1 shows the comparison with 
distances from the surface brightness technique by Gieren 
et al. (2005). Their distances are also larger on average 
than our estimates by 0.15 ±0.09 mag as in the previous 
MS-fitting studies. However, they recalibrated the projec- 
tion factor, which was used to convert observed velocities 
to pulsational velocities, using MS-fitting distances to clus- 
ter Cepheids. Therefore, their estimates are not completely 
independent from the above MS-fitting studies. 

Our short distances are mainly due to our lower E(B-V) 
values than those in the previous work. This is because 
the MS-fitting distances are correlated with reddening by 
A(m-M) /E(B-V) ~ 2 (Paper III). The top panel in Fig- 
ure 12 compares the Cepheid reddening with those in Feast 
& Walker (1987), and individual estimates are listed in Ta- 
ble 1 1 . They estimated the OB star reddening at the location 
of the Cepheids and then transformed it to the value appro- 
priate for the Cepheid using a color-dependent reddening law 
(e.g., eqs. [2] and [3]). The mean difference from our study is 
(E(B-V)} = 0.04±0.02 mag, our values being smaller on av- 
erage. The difference in E(B-V) then approximately matches 
the size expected to produce the difference in distance modu- 
lus (~ 0.1 mag). 

Individual reddening estimates from Fernie et al. (1995), 
Tammann et al. (2003), and Hoyle et al. (2003) are also listed 
in Table 11, and comparisons with our estimates are shown 
in Figure 12. Many studies used an extensive compilation of 
Cepheid reddening in the David Dunlap Observatory (DDO) 
database of Galactic classical Cepheids, 6 which is from het- 
erogeneous methods (Fernie 1990, 1994; Fernie et al. 1995). 
The mean difference with their data is A(E(B -V)) = 0.05 
in the sense that our reddening is smaller. Tammann et al. 
adopted these reddening estimates for their cluster Cepheids 
but made slight corrections to remove a mild correlation of 
the color excess with residuals from the period-color rela- 
tion. This results in a reduced difference with our estimates 
by A(E(B-V)) = 0.02. However, the results from Tammann 
et al. are internally less consistent, in the sense that the extinc- 
tion corrections on the Cepheid magnitudes and the reddening 
values for the MS-fitting distance estimates were not derived 
from the same methods. Finally, we compared our redden- 

6 http://www.astro.utoronto.ca/DDO/research/cepheids/ 
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FIG. 12. — Comparison of Cepheid E(B — V) with those in Feast & Walker 
(1987), Fernie et al. (1995), Tammann et al. (2003), and Hoyle et al. (2003). 
The differences are in the sense of their values minus those in this paper. 
Error bars represent our error estimates in the top three panels, but they show 
a quadrature sum from two studies in the bottom panel (Hoyle et al.). 

ing values with those in Hoyle et al., which were based on 
the cluster UBV color-color diagram. Their values are also 
systematically larger than ours by A(E(B-V)) = 0.04±0.01. 

There are a limited number of reddening estimates from 
the spectral types in the literature. Kraft (1958) estimated 
E(B-V) = 0.45 from one B-type star in NGC 129, while Arp 
etal. (1959) found (E(B-V)) = 0.53 ±0.01 from seven stars 
in the vicinity of its Cepheid (DL Cas). In M25, Sandage 
(1960) estimated (E(B-V)) = 0.43 ±0.05 from the spectral 
types given by Feast (1957) and 0.43 ±0.06 by Wallerstein 
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TABLE 1 1 

Comparison of Cepheid E(B-V) with Published Results 



Cluster 


This Study 11 


FW87 


F95 


TSR03 


HST03 


NGC 129 . 


0.46 ±0.03 


0.49 


0.52 ±0.01 


0.48 


0.52 ±0.05 


NGC 5662 


0.26 ±0.03 


0.29 


0.28 ±0.02 


0.26 




Lynga 6 . . . 


1.23 ±0.03 


1.22 


1.27 ±0.02 


1.21 


1.24 ±0.08 


NGC 6067 


0.32 ±0.03 


0.32 


0.33 ±0.01 


0.32 


0.33 ±0.03 


NGC 6087 


0.12 ±0.03 


0.17 


0.19 ±0.01 


0.18 




M25 


0.39 ±0.03 


0.44 


0.43 ±0.01 


0.40 


0.45 ±0.04 


NGC 7790 b 


0.48 ±0.02 


0.59 


0.59 ±0.02 


0.55 


0.55 ±0.05 



NOTE. — FW87, Feast & Walker (1987); F95, Fernie et al. (1995); TSR03, Tammann et al. (2003); 
HST03, Hoyle et al. (2003). 

"Errors from differential reddening are not included. 
b Average value for CEa Cas, CEb Cas, and CF Cas. 



(1957), but found a slightly larger reddening, 0.46 ± 0.04, in 
the vicinity of its Cepheid (U Sgr). Finally, Kraft (1958) esti- 
mated (E(B-V)) = 0.49 ±0.02 from six stars in NGC 7790. 
These estimates are generally found between the previous es- 
timates and our new values. Because of low Galactic latitudes 
of the clusters (\b\ < 5°), the dust map of Schlegel et al. (1998) 
could not provide meaningful upper limits on the color excess. 

We contend that our reddening estimates are more reliable 
than the previous values, which mostly relied on the solar 
metallicity ZAMS -fitting on UBV color-color diagrams. In 
contrast to these studies, we took into account metallicity and 
age effects on MS fitting, and our estimates are based on si- 
multaneous fits over a wider range of spectral bands, includ- 
ing near-infrared JHK S . We also note that photometric zero- 
point errors in i/-band (e.g., Bessell 2005) could introduce 
significant uncertainties in color excesses derived from UBV 
color-color diagrams. 

5. GALACTIC PERIOD-LUMINOSITY RELATIONS 

With cluster parameters derived in the previous section, 
we are now in a position to estimate absolute magnitudes of 
Cepheids and to redetermine the Galactic P-L relations. We 
show that our new P-L relation in V is generally fainter than 
in previous studies. However, it is in good agreement with 
recent parallax studies (Benedict et al. 2007; van Leeuwen et 
al. 2007) when the Wesenheit magnitudes are employed. 

5.1. Absolute Magnitudes of Cepheids 

The absolute V magnitude was computed for each Cepheid 
from 

(M v ) = (m v )-(m-M)o-R v E(B-V), (10) 

where (m v ) is the observed intensity-mean magnitude and 
(m-M)o is the extinction-corrected distance modulus. We 
used the same color-dependent prescriptions for Ry and E(B- 
V) as those for MS dwarfs (eqs. [2] and [3]). The absolute 
magnitudes in other filter passbands were computed from the 
CCM89 reddening law at each cluster's R V fl- Table 12 lists ab- 
solute magnitudes in BVI C JHK S for all Cepheids in this study. 
In the last column we also show the reddening-insensitive We- 
senheit magnitude, (W VI ) = (V)-2A5((V)-(I C )), which was 
adopted in the HST Key Project. 7 

7 We simply adopted A v /E(V-I C ) = 2.45 from Freedman et al. (2001) to 
directly compare our results with those of other studies, which also employed 
the same relation. We would derive Ay /E{V -I c ) = 2.56 for the Cepheids 



Table 13 shows systematic errors in (My) of the Cepheids. 
As in Table 8, we estimated the size of errors by perform- 
ing MS fits for alternate values of systematic error quantities. 
These are shown in the first and second columns, which are 
the same as those in Table 8 except the last two rows: Amy 
and <r E (B-v)- The Amy represents a systematic error in the 
Cepheid photometry, which was treated as independent of the 
cluster photometry and any of the MS-fitting process. Note 
that (m-M)o, Rv, and E(B-V) in equation (10) are corre- 
lated with each other via MS fitting. 

The (Je(b-v) in Table 13 represents the error in the redden- 
ing for individual Cepheids. In equation (10) we adopted the 
Cepheid reddening as the mean cluster reddening value. How- 
ever, reddening for individual stars becomes more uncertain 
as we deal with more differentially reddened clusters. For 
these clusters, the colors for individual stars scatter along a 
reddening vector, and this can be detected by a large size of 
the MS width above photometric precision. 

To estimate the size of differential reddening for each clus- 
ter, we first computed the standard deviation of B - V colors 
in the upper MS from Table 7. Since the photometric filter- 
ing likely reduces the size of the dispersion, we corrected for 
it from artificial cluster CMD tests (Paper III). We generated 
solar metallicity, 80 Myr isochrones with photometric errors 
of 0.02 mag in colors and magnitudes. We used the Salpeter 
(1955) mass function for the primaries and a flat mass func- 
tion for the secondaries for a 40% binary fraction. 8 Single 
stars and binaries were then randomly displaced from the MS 
assuming a normal distribution of individual reddening with 
a standard deviation <t E (b-v)- After applying the photometric 
filtering, we then estimated the size of the reduction in disper- 
sion as a function of oe(b-v)- In this way, we inferred oe{b-v) 
for each cluster from the observed dispersion in B - V , which 
are shown in the second column of Table 14. We made no 
corrections to those for NGC 5662 and NGC 7790 because 
the dispersions in B - V were equal to or smaller than our as- 
sumed photometric precision (<t B -v = 0.02). 

In addition to previously known clusters with differential 
reddening (NGC 129, NGC 5662, and M25), we also found 

with (Bo) _ (^o) = 0.7 at Rv = 3.35, computed from the average Ry,o of our 
cluster sample. This difference in Ay /E(V - Ic) has a negligible effect on 
Wvi of Cepheids. It is noted that Macri et al. (2001b) tested the CCM89 
reddening law using VIH photometry for 70 extragalactic Cepheids in 13 
galaxies. They showed that a mean color-excess ratio of E(V -H)/E(V -I) 
is in good agreement with the CCM89 predicted value. 

8 The binary fraction is defined as the number of binaries divided by the 
total number of systems in the considered fitting range. 
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TABLE 12 
Cepheid Absolute Magnitudes 



Cepheids 


logP 


(M B ) 


(M v ) 


(Mi c ) 


(Mj) 


(M„) 




(M W (vi)) a 


DL Cas . . 


0.903 


-2.988 


-3.673 


-4.364 








-5.285 


V Cen ... 


0.740 


-2.754 


-3.371 


-4.042 


-4.557 


-4.848 


-4.951 


-4.974 


TWNor . 


1.033 


-3.187 


-3.957 


-4.716 


-5.287 


-5.580 


-5.775 


-5.676 


V340 Nor 


1.053 


-2.879 


-3.720 


-4.505 


-5.142 


-5.500 


-5.641 


-5.604 


SNor ... 


0.989 


-2.864 


-3.685 


-4.518 


-5.128 


-5.470 


-5.593 


-5.688 


USgr ... 


0.829 


-2.829 


-3.535 


-4.270 


-4.794 


-5.086 


-5.196 


-5.294 


CEaCas . 


0.711 


-2.458 


-3.129 


-3.951 








-5.090 


CEbCas . 


0.651 


-2.306 


-2.998 


-3.731 








-4.739 


CFCas .. 


0.688 


-2.189 


-2.911 


-3.667 








-4.707 



"Wesenheit index from Freedman et al. (2001): (W(VI)) = (V) -2.45((V) - (/c». 



TABLE 13 
Error Budget in (M v ) 



Source of Error 


AQuantity 


NGC 129 


NGC 5662 


Lynga 6 


NGC 6067 


NGC 6087 


M25 


NGC 7790 




DL Cas 


VCen 


TWNor 


V340 Nor 


SNor 


USgr 


CEa, CEb, CF Cas 


|\1 




^0.032 


=f0.028 


=F=0.056 


=1=0.054 


=1=0.028 


T-0.045 


=1=0.054 


log; (Myr) 


±0.2 


±0.037 


±0.007 


=F=0.014 


±0.057 


±0.061 


±0.051 


±0.021 


Helium (AY) . . . 


±0.010 


T-0.030 


T0.030 


=1=0.030 


=1=0.030 


=1=0.030 


=1=0.030 


=1=0.030 


A(B-V) 


±0.025 


T0.009 


=fo.ooi 


±0.014 


=1=0.007 


=1=0.002 


=1=0.002 


±0.001 


A(V-I) C 


±0.025 








=1=0.012 






=p0.011 


Al 


±0.025 


±0.048 


±0.025 


±0.028 


±0.030 


±0.022 


±0.013 


±0.026 


AJ 


±0.011 


=1=0.014 


=1=0.004 


=f0.002 


=1=0.003 


=1=0.005 


=1=0.005 


±0.000 


AH 


±0.007 


±0.001 


=1=0.001 


=F=0.002 


=1=0.001 


=1=0.002 


=1=0.002 


±0.001 


A A'. 


±0.007 


±0.010 


±0.002 


=F=0.002 


±0.002 


±0.004 


±0.004 


=1=0.001 


Fitting 




±0.023 


±0.044 


±0.094 


±0.032 


±0.031 


±0.040 


±0.023 


Am v b 


±0.025 


±0.025 


±0.025 


±0.025 


±0.025 


±0.025 


±0.025 


±0.025 


VE(B-V) C 




±0.347 


±0.066 


±0.316 


±0.159 


±0.197 


±0.379 


±0.041 


Total 




±0.357 


±0.097 


±0.338 


±0.187 


±0.216 


±0.390 


±0.089 



NOTE. — A(B— V), A(V-T)c, AV, AJ, AH, and AK S represent zero-point errors in the cluster photometry. 
"Individual errors taken from Table 1 . 

b Photometric zero-point error in the intensity mean magnitude of a Cepheid. 
c Differential reddening over a cluster field. 



dispersion in Table 14 as a characteristic measure of differen- 
tial reddening as a conservative error estimate. The <Je(b-v) 
in Table 14 was then multiplied by R v to estimate the error 
in (My), which is shown in Table 13. For all clusters except 
NGC 7790, differential reddening is the largest error source 
in (Mv). 

In previous studies, local Cepheid reddening was often de- 
termined using photometry of neighboring stars. We exper- 
imented with this approach for M25. We estimated individ- 
ual Cepheid's reddening on (B — V,V) and assumed that the 
dispersion around the best-fitting isochrone is solely due to 
differential reddening. However, we could not place a strong 
constraint on the local color excess for U Sgr because there 
were not a sufficient number of stars with good E(B-V). In 
addition, cluster binaries cannot be easily distinguished from 
highly reddened single stars on CMDs, which could lead to 
an overestimation of reddening. 

5.2. The Galactic P-L Relations 

Our Cepheid sample does not span a wide range of periods, 
which leads to a large error in P-L slopes. Instead, we adopted 
the P-L slopes derived from LMC Cepheids, assuming that 
the P-L relations for Cepheids in the Galaxy and in the LMC 
have the same slopes. Some investigators have claimed that 



TABLE 14 
Results from Artificial Cluster 
Tests 



Cluster 


<?E(B-V) d 


A{M v ) h 


NGC 129 . 


0.105 


+0.04 ±0.14 


NGC 5662 


0.020° 


-0.06 ±0.05 


Lynga 6 . . . 


0.100 


+0.04 ±0.14 


NGC 6067 


0.050 


-0.06 ±0.07 


NGC 6087 


0.055 


-0.06 ±0.08 


M25 


0.120 


+0.02 ±0.15 


NGC 7790 


0.013 c 


-0.05 ±0.04 



a Size of differential reddening inferred from arti- 
ficial cluster tests, assuming a normal distribution of 
E(B-V). 

b This value should be subtracted from the esti- 
mated (M v ). Errors are semi-interquartile ranges. 
c Not corrected from the artificial cluster tests. 



that the remaining clusters except NGC 7790 have larger dis- 
persions than photometric errors. While this could be due to 
remaining binaries, background stars, rapid rotators, or an un- 
derestimation of the photometric errors, we took the excess 
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the Galactic P-L relations are steeper than those of the LMC 
Cepheids, and that the LMC P-L relations have a break at P w 
10 days, possibly due to different metal contents of Cepheids 
(Kanbur & Ngeow 2004; Sandage et al. 2004; Ngeow et al. 
2005). However, Macri et al. (2006) recently showed that P- 
L relations from the LMC Cepheids in the OGLE-II catalog 
(Udalski et al. 1999a) are a good fit to Cepheids in the two 
fields of NGC 4258, which have approximately the LMC and 
the Galactic metal abundances (see § 6.1). 

We derived LMC P-L slopes using fundamental mode 
Cepheids in the OGLE-II catalog (Udalski et al. 1999a) and 
those in Persson et al. (2004). The OGLE-II Cepheids have an 
average period of (logP(days)) w 0.6 with about half of them 
being in the same period range as our Galactic Cepheid sam- 
ple. The P-L relations in the LMC from the OGLE-II database 
were originally derived by Udalski et al. (1999b) and then re- 
vised according to the new photometric calibration (Udalski 
2000). We rederived the period-apparent magnitude relations 
in BVIc using the revised data from the OGLE-II Internet 
archive 9 after iterative 2.5cr rejection from 680 Cepheids with 
logP > 0.4 (Udalski et al. 1999b). We also derived similar 
relations in JHK S on the LCO system from the photometry in 
Persson et al. (2004). As in their analysis, we excluded four 
stars (HV 883, HV 2447, and HV 2883 have periods longer 
than 100 days, and HV 12765 is about 0.2 mag brighter than 
the mean P-L relation), leaving 88 Cepheids with logP > 0.4. 

We then derived the Galactic P-L relations in BVIcJHK, 
after correcting for interstellar extinction: 

(M B ) = (-3.169±0.085)-(2.503± 0.048)(logP-l), 

(M v ) = (-3.932±0.066)-(2.819±0.032)(logP-l), 

(M /c ) = (-4.712±0.057)-(3.004±0.021)(logP-l), 

(My) = (-5.271 ±0.076)-(3. 148 ±0.053)(logP- 1), 

(M H ) = (-5.593 ±0.069)-(3.233±0.044)(logP-l), 

(M Ks ) =(-5.718 ±0.064)- (3.282 ±0.040)(logP- 1). 

The slope errors are those estimated from the LMC Cepheids. 
The zero-point errors were estimated from the magnitude dis- 
persion among the Cepheids, assuming that the magnitude er- 
rors are uncorrected with each other. We took the average 
magnitudes and periods of the three Cepheids in NGC 7790 
as one data point. Cepheids in the two northern clusters 
(NGC 129 and NGC 7790) have no photometry in Laney & 
Stobie (1992), so the P-L relations in JHK S were derived from 
five clusters. 

Figures 13 and 14 display P-L relations for our Galactic 
Cepheids {bull's-eyes) in all filter passbands. In V-band the 
X 2 of the fit is 8.0 for 6 degrees of freedom, which indicates 
our reasonable error estimation. The standard deviation of 
the Galactic Cepheids around the P-L relation is 0.18 mag, 
which is comparable to that of the LMC counterparts (Udalski 
et al. 1999b). This value is smaller than previous estimates 
for the Galactic Cepheids from MS fitting: Tammann et al. 
(2003) found 0.26 mag for 25 Cepheids (see also Sandage et 
al. 2004), and Fouque et al. (2003) estimated 0.27 mag from 
24 Cepheids. 

In Figure 13 the OGLE-II Cepheids are also shown at a dis- 
tance modulus of 18.50 mag, which was adopted in the HST 
Key Project (small boxes). Magnitudes of these Cepheids 
were corrected for extinction based on the OGLE-II reddening 
map, which was derived from the / c -band magnitudes of red 
clump stars (Udalski et al. 1999a). As seen in the figure, the 

9 http://w w w. astrouw. edu. pl/~ogle/ogle2/cep_lmc .html 
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FIG. 13. — Galactic P-L relation in V. The bull's-eyes are our Galactic cali- 
brators with their lcr errors, and the solid line is a fit to these points assuming 
the LMC P-L slope. The small boxes are LMC fundamental mode Cepheids 
from the OGLE-II survey, assuming (m-M) = 18.50 and the OGLE-II red- 
dening map. The dashed line is a fit to these points. 

LMC P-L relation (dashed line) is brighter than our Galactic 
P-L relation (solid line) by <~ 0.3 mag. This is due to metal- 
licity and reddening effects, as well as the LMC distance as 
discussed in § 6.2. 

5.3. Comparison with Previous Studies 

Figure 15 compares our (My) with those of Feast & Walker 
(1987), Tammann et al. (2003), Hoyle et al. (2003), and 
Gieren et al. (2005). Error bars represent our error estimates 
only. We revised absolute magnitudes in Feast & Walker and 
Tammann et al. assuming our Pleiades distance (5.63 mag 
in distance modulus). Individual (My) estimates are shown 
in Table 15. We restricted our comparison to (My) since 
absolute magnitudes in other passbands are correlated with 
each other. The weighted average differences are A(My) = 
-0.17 ±0.07, -0.09 ±0.09, -0.14 ±0.10, and -0.22 ±0.11 
from Feast & Walker, Tammann et al., Hoyle et al., and Gieren 
et al., respectively, in the sense that our estimates are fainter 
on average. This is mainly due to the fact that our distance 
and reddening estimates are smaller on average than those in 
the previous studies (§ 4.2). 

In addition, we also made a comparison with all cluster 
Cepheids in Tammann et al. (2003) as shown in the left panel 
of Figure 16. As found in the above Cepheid-to-Cepheid com- 
parison, the difference with their (My) is statistically insignif- 
icant for Cepheids with logP < 1 . 1 . However, their P-L rela- 
tion is steeper than ours because their long-period Cepheids at 
logP > 1.1 are brighter than expected from our P-L relation. 
If we assume their P-L slope for our sample, we would derive 
(M v ) = -4.01 for 10-day period Cepheids, which is 0.25 mag 
fainter than their P-L zero point. 

In the right panel of Figure 16, we compare our Galac- 
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FIG. 14.— Same as Fig. 13, but in all BVI C JHK S . 



TABLE 15 

Comparison of Cepheid (M v ) with Published Results 



Cluster This Study FW87 a TSR03 a HST03 G05 



DLCas -3.67 ±0.36 -3.92 -3.77 -3.69 

VCen -3.37±0.10 -3.30 -3.19 ••■ -3.45 

TW Nor -3.96 ±0.34 -3.85 -3.72 -3.79 

V34()Nor -3.72±0.19 -3.86 -3.85 -3.92 -3.91 

SNor -3.69 ±0.22 -4.01 -4.00 ••■ -4.21 

USgr -3.54±0.39 -3.83 -3.66 -3.88 -3.62 

CEa,CEb,CF Cas b -3.01 ±0.09 -3.36 -3.37 -3.34 



NOTE. — FW87, Feast & Walker (1987); TSR03, Tammann et al. (2003); HST03, Hoyle et al. 
(2003); G05, Gieren et al. (2005). 

"Magnitudes revised assuming the same Pleiades distance, (m-M)o = 5.63, as in this paper. 
b Average value for CEa Cas, CEb Cas, and CF Cas. 



tic P-L relation in V with that of the recent HST parallax 
study (Benedict et al. 2007), which provided accurate paral- 
laxes (average 8%) for nine Galactic Cepheids using the HST 
Fine Guidance Sensor (see also Benedict et al. 2002b). Their 
Cepheid parallaxes were estimated with respect to photomet- 
ric parallaxes of reference stars, which set the absolute dis- 
tance scale of each reference frame. The magnitudes were 
then corrected for interstellar extinction from either color- 
color relations or the DDO Galactic Cepheid database. We 
have no Cepheids in common with their study since their dis- 
tance measurements were focused on nearby field Cepheids. 
If we assume their P-L slope, our 10-day period Cepheids 
would be 0.20 mag fainter than their zero point. If the differ- 
ence is solely due to the error in our reddening estimates, our 
estimates would have been underestimated by A(E(B— V)} <~ 
0.04. However, it is noted that reddening can be determined 
far more accurately for cluster Cepheids than it can for field 
stars. Our error bars in Figure 16 are slightly larger than those 



of the HST study because of our conservative error estimates 
for the extinction correction. 

On the other hand, the ffi>r-based P-L relation is in better 
agreement with ours when the Wesenheit magnitude W(VI) is 
used as shown in Figure 17. The weighted mean difference 
of our Galactic Cepheids from their best-fitting P-L relation 
(solid line) is A(Mu/(v/)) = 0.05 ±0.06 in the sense that our 
P-L relation is fainter. Since the Wesenheit magnitude was 
designed to minimize the effect of errors in the reddening, the 
agreement may indicate that the large difference in (My) is 
due to the difference in the distance estimation. 

Previous studies based on Hipparcos data typically yielded 
(My) Ri -4.2 for 10-day period Cepheids (Feast & Catchpole 
1997; Feast et al. 1998; Lanoix et al. 1999; Groenewegen & 
Oudmaijer 2000), while Luri et al. (1998) obtained (M v ) w 
-3.9 from statistical parallaxes. Recently, van Leeuwen et 
al. (2007) reexamined the Galactic P-L relations based on 
the new reduction of the Hipparcos astrometric data (van 
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FIG. 15. — Comparison of Cepheid absolute magnitudes with those in Feast 
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TABLE 16 

Error Budget in (M v ) for 10-day Period 
Cepheids from MS Fitting 




Source of Error 


A{M V ) 


Isochrone Calibration 


±0.02 


Metallicity Scale (±0.1 dex) 


±0.11 


Age Scale [Alog?(Myr) = ±0.2] . . 


±0.03 


Reddening Laws 


±0.05 


Differential Reddening + Binary 


±0.05 


Total Systematic 


±0.14 


Statistical 


±0.07 


Combined (Systematic + Statistical) 


±0.16 



Leeuwen & Fantino 2005), where improvements of up to a 
factor of 2 in parallax accuracy have been achieved by recon- 
structing the satellite's attitude (van Leeuwen 2005). They 
analyzed their data together with the HST parallaxes to de- 
rive a P-L relation in W(VI). If we adopt their P-L slope in 
W(VI), our P-L relation would be A(M wm ) = 0.07 ±0.06 
fainter than their best-fitting P-L relation. 

5.4. Systematic Errors in P-L relations 

When determining the P-L relations in the previous section, 
we assumed that the errors in (My) of individual Cepheids 



(Table 13) are uncorrected with each other. However, there 
are a number of error sources that could systematically change 
(My). These are shown in Table 16 with their estimated size 
of error contributions to the P-L zero point, which is described 
below in detail. In short, the calibration error is from the un- 
certainty in the adopted distance of the Pleiades (Paper III). 
Errors in the metallicity scale of ±0. 1 dex and the age scale of 
Alogf(Myr) = 0.2 were adopted. The error from the redden- 
ing laws represents the case of using different color-dependent 
prescriptions for E(B-V) andR v from Feast & Walker (1987) 
and Laney & Stobie (1993). We performed artificial cluster 
CMD tests to estimate a probable size of a bias in the MS- 
fitting technique from differential reddening and unresolved 
cluster binaries. 

From the comparison with Yong et al. (2006) and An- 
drievsky et al. (2002), we inferred that our adopted metallic- 
ity scale (Fry & Carney 1997) has a 0.1 dex systematic error 
(§ 2.1). Metallicity changes of ±0.1 dex would vary the (My) 
of 10-day period Cepheids from -3.84 to -4.04 mag. 

In this paper, only the relative age scale is relevant with re- 
spect to the Pleiades. Our assumed age range in § 2.1 reflects 
the cluster-to-cluster variations with respect to the adopted 
age of 80 Myr. However, it is possible that the Cepheid clus- 
ters are all systematically younger or older than the Pleiades. 
This would change the reddening values since the mean color 
of the upper MS varies with the cluster age. In fact, some 
of the CMDs have stars that extend far above the top end of 
the 80 Myr isochrone. If these bright stars are cluster mem- 
bers, they would suggest younger ages than 80 Myr. Never- 
theless, the P-L zero point is relatively insensitive to the age. 
If the overall age scale is uncertain at the la level of individ- 
ual cluster ages (50-120 Myr), we would have (My) = -3.90 
and -3 .96 mag for 1 0-day period Cepheids at 50 and 1 20 Myr, 
respectively. 

We estimated the error from the reddening laws using dif- 
ferent prescriptions for the color-dependent E(B-V) and Ry 
in Feast & Walker (1987) and Laney & Stobie (1993). For the 
same E(B - V)o and Ry.o values derived from the cluster MS 
fitting, we found that these prescriptions result in a fainter P-L 
relation by A(M V ) w 0.05 than our default case (eqs. [2] and 
[3]). 

The average Ry.o that we inferred from the clusters in this 
study is 3.19 (Table 9), which is in good agreement with the 
average Galactic value for the diffuse interstellar medium, 
Rv ~ 3.1. The formal standard deviation of Ry,o is 0.17, but 
individual errors are too large to detect a cluster-to-cluster 
variation. Even if we adopt a different Ryfl, it would have 
a reduced impact on (My). This is because a larger Ry makes 
an MS-fitting distance shorter (Paper III), but it makes the ex- 
tinction correction larger at the same time (eq. [10]). 

In § 5.1 we considered the case where the presence of dif- 
ferential reddening makes it difficult to estimate individual ex- 
tinctions for Cepheids. This is independent of the MS-fitting 
technique in the sense that it does not affect any of the cluster 
parameters. However, differential reddening could also mod- 
ify cluster parameters because of cluster binaries. Unresolved 
binaries are typically brighter and redder than a single-star 
MS, which would make a distance apparently shorter. If there 
exists differential reddening across the cluster field, the distri- 
bution of these binaries on CMDs would be modified, and the 
MS-fitting parameters would subsequently change. 

To estimate the effect of differential reddening, we per- 
formed artificial cluster CMD tests as described earlier in 
§ 5.1. We assumed that E(B-V) has a normal distribution 
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with a standard deviation <j E{B _ v) . For each set of artificial 
cluster CMDs we applied the photometric filtering and esti- 
mated R V fl, E(B-V)o, (m-M)o, and (M v ) as for the actual 
cluster data, assuming (Bo) - (Vq) = 0.7 for Cepheids. Each 



panel in Figure 18 shows the bias in each of these quantities 
as a function of <j E (b-v)- The thick solid line shows the median 
of the bias from 200 artificial cluster CMDs, and each set of 
CMDs contained 60 single stars and 40 unresolved binaries 
in our MS-fitting range for the distance determination. Thin 
lines on either side are the first and third quartiles of these 
estimates. The dashed line indicates a zero bias. 

Even in the presence of a small differential reddening, low 
mass-ratio binaries cannot be easily detected from the photo- 
metric filtering. Because these unresolved binaries are red- 
der than a single-star MS at a given M v , the reddening would 
be overestimated as shown in Figure 18 for (Je(b-v) J$ 0.05. 
A higher E(B - V) then results in a longer distance, and a 
Cepheid would look brighter by several hundredths of a mag- 
nitude in (My). As differential reddening becomes stronger, 
binaries would be more scattered around the MS, and it would 
become further difficult to distinguish them from highly red- 
dened single stars. As a result, both cluster reddening and 
distance estimates decrease at <Je(b-v) ^ 0.05, while semi- 
interquartile ranges continuously increase. We also exper- 
imented with 30% and 50% binary fractions but found in- 
significant changes within the semi-interquartile ranges. 

The last column in Table 14 shows the size of the bias in 
(My) for each Cepheid, which was determined at oe{b-v) in 
the second column. Applying these corrections to (My) in 
Table 12 leads to 0.05 mag reduction in (M v ) for 10-day pe- 
riod Cepheids, again making worse the agreement with the 
brighter calibration in the previous studies. The total % 2 of 
the fit becomes smaller (6.9 for 6 degrees of freedom) if we 
add in quadrature the semi-interquartile range to the error in 
{M v ). 

In summary, our Galactic P-L relation in V-band has a sta- 
tistical error of 0.07 mag and a systematic error of 0. 14 mag in 
the zero point (Table 16). If we add these quantities in quadra- 
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ture, the combined error would be 0.16 mag. Since we deter- 
mined the Galactic P-L relations in other filter passbands with 
the same cluster distance, reddening, and R v estimates as in V, 
their zero-point errors are correlated with each other. There- 
fore, we took 0.16 mag as a systematic error in the Cepheid 
distance scale in the following section. 

6. EXTRAGALACTIC DISTANCE SCALES 

Cepheids have long served as our most commonly used 
standard candles for extragalactic distance studies. The con- 
ventional way of estimating extragalactic Cepheid distances is 
to use a fiducial P-L relation derived from the LMC Cepheids 
with an adopted distance to the LMC. However, geometric 
distances to a few galaxies have become available in recent 
years, providing an opportunity to compare with the Cepheid 
distance scale. In this section we therefore apply our Galactic 
P-L relations to estimate a distance to NGC 4258 (Ml 06) that 
has an accurate geometric distance measurement from its wa- 
ter maser sources (Herrnstein et al. 1999). We then infer a dis- 
tance to the LMC and compare our result with those adopted 
by the HST Key Project and the SN la calibration program 
(Sandage et al. 2006), from which we consider a possible in- 
crease in their Hubble constant Ho measurements. We finish 
with a discussion of the potentially interesting case of M33 
(NGC 598). Our most fundamental conclusion is that the dis- 



TABLE 17 
Distance Modulus and Metallicity 
Sensitivity in NGC4258 



Filter 


(m- 


-M) 


7 

mag dex~' 


Inner Field 


Outer Field 


B 


29.25 ±0.05 


29.33 ±0.07 


-0.20 


V 


29.26 ±0.04 


29.32 ±0.05 


-0.16 


k .... 


29.26 ±0.03 


29.32 ±0.04 


-0.17 


Average 


29.26 ±0.02 


29.32 ±0.03 


-0.17 



NOTE. — Distances are shown at the best-fit reddening from 
the Galactic P-L relations (see text). Errors represent fitting un- 
certainties, and do not include the error from the correlation with 
reddening. 



tance to NGC 4258 is in good agreement with the geometric 
measurement when using the P-L relations inferred from open 
clusters. We also compare our LMC distance with those from 
the HST (Benedict et al. 2007) and the Hipparcos parallax 
studies (van Leeuwen et al. 2007). 

6.1 . The Maser-host Galaxy NGC 4258 

Herrnstein et al. (1999) inferred a geometric distance to 
NGC 4258 from the orbital motions of water maser sources on 
its nucleus and found (m-M) = 29.29 ± 0.09 (statistical) ± 
0.12 (systematic) (see also Argon et al. 2007). Many 
Cepheids were also observed in this galaxy, providing an op- 
portunity to check the Cepheid distance scale (Maoz et al. 
1999; Newman et al. 2001). Recently, Macri et al. (2006) dis- 
covered 28 1 Cepheids and provided accurate BVI photometry 
using the Advanced Camera for Surveys/Wide Field Camera 
(ACS/WFC) onboard the HST. In particular, they observed 
two fields in the galaxy, located at two different galactocentric 
distances with significantly different gas-phase metal abun- 
dances, and derived the metallicity sensitivity of the Cepheid 
luminosity. 

We reexamined the Macri et al. data set in light of our 
revised P-L relationship. Following their selection proce- 
dure, we used the "restricted" sample of 69 Cepheids and es- 
timated distance moduli for the two groups of Cepheids in 
the inner (with a period cut of P > 12 days) and outer fields 
(P > 6 days). For each group of this sample we corrected 
for extinction using our P-L relations in BVIc- The abso- 
lute distance to NGC 4258 was then estimated by anchoring 
our Galactic sample to a reference gas-phase metal abundance 
within the galaxy. 

Figure 19 shows how we determined the average distance 
and reddening. In the left two panels apparent distance mod- 
uli in BVIc are shown for Cepheids in the inner and outer 
fields {top and bottom, respectively). Distances from shorter 
wavelength filters are systematically longer, and the data are 
well fitted by a CCM89 extinction curve (solid line). Here we 
assumed Ry = 3.35 for the Cepheids, which was computed 
from the average R v ,o of our cluster sample with equation (2) 
at (B) -{V)o = 0.72. At this value, the CCM89 extinction 
law yields absorption ratios A\/Ay of 1.30 and 0.61 forB and 
Ic, respectively, where we took effective wavelengths of these 
filters from BCP98. We added 0.02 mag in quadrature to the 
fitting error to take into account the photometric zero-point 
error. At the limit of A" 1 — > 0, an intercept on the ordinate 
is the true or unreddened distance modulus (e.g., Madore & 
Freedman 1991, 1998). The right two panels in Figure 19 
show likelihood contours in the average reddening and dis- 
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FIG. 19. — Top left: Apparent distance moduli of inner field Cepheids in NGC 4258 with the best-fitting CCM89 extinction curve (solid line). Top right: 
Likelihood contours in the average reddening and distance modulus of inner field Cepheids shown at A\ = 2.30, 6.17, and 11.8 (68.3%, 95.4%, and 99.73% 
confidence levels for 2 degrees of freedom). The horizontal lines represent a geometric distance from the water maser sources and its \o error (Herrnstein et al. 
1999). Bottom left: Same as the top left panel, but for outer field Cepheids. Bottom right: Same as the top right panel, but for outer field Cepheids. 



tance modulus from these fits. Contours are shown at A% 2 = 
2.30, 6.17, and 11.8 (68.3%, 95.4%, and 99.73% confidence 
levels) for 2 degrees of freedom. Horizontal lines represent 
the distance estimate from the water maser sources and its 
combined lcr error (Herrnstein et al. 1999). Distance moduli 
at the best-fit reddening are listed in Table 17 for each of the 
two galactic fields. 

When the reddening is constrained in this way, distances 
in both fields are in good agreement with the maser distance. 
However, the distance from the inner field is shorter than the 
distance from the outer field by A(m-M)o = 0.07. The most 
natural interpretation of this difference is that it measures the 
metallicity dependence of the Cepheid luminosity. This effect 
has been extensively discussed in the literature and tested em- 
pirically (Freedman & Madore 1990; Gould 1994; Sasselov 
et al. 1997; Kochanek 1997; Kennicutt et al. 1998; Udalski 
2000; Kennicutt et al. 2003; Groenewegen et al. 2004; Sakai 
et al. 2004; Storm et al. 2004; Romaniello et al. 2005; Macri 
et al. 2006). 

To estimate the true distance modulus of the galaxy, we 
first estimated metal abundances in the two NGC 4258 fields 
and then inferred a distance at the metallicity of the Galactic 
Cepheids. Stellar abundances cannot be measured directly for 
most of the Cepheid galaxies. Instead, the HST Key Project 
adopted "empirical" gas-phase [O/H] abundances from Zarit- 
sky et al. (1994). Macri et al. computed oxygen abundances at 
each Cepheid's location from the H II region abundance gradi- 
ent in Zaritsky et al.. The average abundances from this table 
are (12 + log (O/H)) = 8.94 and 8.54 dex for the inner and 
outer field Cepheids, respectively, when the LMC Cepheids 



have 12 + log (O/H) = 8.50 on this scale (Sakai et al. 2004). 
To relate these abundances to those of our Galactic Cepheid 
sample, we adopted a stellar abundance difference between 
the LMC Cepheids and Galactic supergiants from Luck et 
al. (1998) of A ([Fe/H]) (Galaxy -LMC) w 0.30. They also 
found a similar a-element enhancement between the LMC 
and the Galactic samples. In the following discussion we 
adopted A (log (O/H)) (Galaxy -LMC) = 0.30 ±0.06, assum- 
ing that the stellar abundance difference is the same as the 
difference in the gas-phase oxygen abundance. 

Since the metallicity of the Galactic Cepheids is found be- 
tween the metal abundances for the inner and outer fields, the 
true distance would lie between distance estimates from the 
two galactic fields. From a linear interpolation, we obtained 
(m-M)o = 29.28 ± 0.10 ± 0. 16 (P-L zero point) for a distance 
to NGC 4258, which is in good agreement with the current 
maser distance. The first error is from the fit to the HST ACS 
data (0.10 mag), an error in R v (+~0.03 mag for ±0.5 change 
in Ry), and an error in the metallicity scale (+~0.02 mag for 
±0.10 dex). We adopted an error in the metallicity scale of 
±0.10 dex to account for a possible difference in the stellar 
and gas-phase abundances, but it has only a moderate effect 
on the distance determination. 

Previous studies found somewhat longer Cepheid dis- 
tances to NGC 4258 although these values are consistent 
within the errors with the maser solution. As a part of the 
HST Key Project, Maoz et al. (1999) estimated (m-M) = 
29.54 ±0.12 ±0.18 using the Wide-Field Planetary Cam- 
era 2 (WFPC2), and Newman et al. (2001) revised it to be 
(m-M)o = 29.47 ±0.09 ±0.15, where the first errors repre- 
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sent those unique to each determination and the second errors 
are systematic uncertainties in the HST Key Project distances. 

However, the reddening in the HST Key Project is signifi- 
cantly larger than our estimates. The final result in the HST 
Key Project lists (E(V-I) C ) = 0.22 ±0.04 or (E(B-V)) = 
0. 17 ± 0.03 for Cepheids located in the middle of galactocen- 
tric distances to the two ACS fields. On the other hand, we 
found (E(B - V)} = 0.08 ± 0.04 and 0.06 ± 0.05 for the inner- 
and outer-field Cepheids, respectively, where the errors are 
the lcr contour sizes in Figure 19. Macri et al. also estimated 
individual Cepheid reddenings using the OGLE-II LMC P-L 
relations, and their average values are also larger than ours, 
being (E(B-V)) = 0.18 and 0.13 mag for the inner and outer 
fields, respectively. 

In Table 17 we also listed the metallicity sensitivity (7) in 
each band, which was computed from the distances and the 
average abundance difference in the two galactic fields. The 
magnitude of this effect is smaller than the original estimate 
in Macri et al., 7 = -0.29 ± 0.09 (random) ± 0.05 (systematic) 
mag dex" 1 . Since our estimate neglects the change in color 
with metallicity, we employed the metallicity corrections 
from Kochanek (1997) in the next two sections. Kochanek 
derived corrections in UBVRIJHK based on a simultaneous 
analysis of Cepheids in 17 galaxies with a proper treatment of 
correlations between extinction, temperature, and metallicity. 
The corrections are +0.20 ± 0. 1 8, -0.08 ± 0. 14, -0.21 ± 0. 14, 
-0.26 ±0.13, -0.34 ±0.14, and -0.40 ±0.13 mag dex" 1 in 
BVIcJHK, respectively, in the sense that metal-rich Cepheids 
have a decreased flux in B but increased fluxes in VIcJHK, 
possibly due to line-blanketing and back-warming among 
other effects. These estimates are based on (E(B-V)} = 0.15 
for LMC Cepheids, but their dependence on the LMC red- 
dening is small. If we adopt these corrections for NGC 4258, 
we would derive (m-M)o = 29.34 ±0.1 3 for the inner-field 
and 29. 17 ± 0. 19 mag for the outer-field Cepheids, where the 
errors represent the la statistical uncertainty. 

6.2. LMC 

Our Galactic P-L relations can also be applied to estimate 
the distance to the LMC and its mean reddening. We included 
JHK S photometry of the LMC Cepheids from Persson et al. 
(2004), as well as BVI C photometry from the OGLE-II cata- 
log in the multi-wavelength fitting technique. This approach 
assumes that Cepheids in these two surveys were drawn from 
the same population. 10 

Figure 20 shows likelihood contours in the solution for 
the average LMC reddening and its distance modulus (mid- 
dle panel) from the fits on the apparent distance moduli in 
BVIcJHK, (left panel). We applied metallicity corrections 
for A (log (O/H)) (Galaxy- LMC) = 0.30 and assumed that 
the metallicity sensitivities on the Glass-Carter JHK system 
in Kochanek (1997) are the same as in the LCO JHK S sys- 
tem. As in the NGC 4258 case, we assumed R v = 3.35 for 
the Cepheids, which yields absorption ratios A\/A v of 0.29, 
0.18, and 0.12 for JHK S , respectively, from the CCM89 ex- 
tinction law. In this calculation we took the effective filter 
wavelengths from Persson et al. (2004) for JHK S on the LCO 
system. We added 0.02 mag in BVI C (Udalski et al. 1999a) 
and 0.03 mag in JHK S (Persson et al. 2004) in quadrature to 
the fitting errors to account for photometric zero-point errors. 

10 The OGLE-II Cepheids were observed along the LMC bar, while the 
Cepheids observed in JHK S (Persson et al. 2004) are more dispersed over the 
LMC. 



TABLE 18 
Distance Modulus of the LMC 



Filter (m-M)o 

OGLE-II E(B- VY Best-fit E(B-Vf 



B 18.19±0.06 18.35±0.06 

V 18.23 ±0.05 18.34 ±0.05 

I c .... 18.29 ±0.05 18.35 ±0.05 

/ ••• 18.35 ±0.05 

H ••• 18.32 ±0.05 

K, .... ■■■ 18.36 ±0.05 

Average 18.24 ±0.03 18.34 ±0.02 



NOTE. — After applying metallicity corrections from 



Kochanek (1997). 

"Distances derived assuming the OGLE-II reddening map. 

b Distances shown at the best-fit reddening from the Galactic P- 
L relations (see text). Errors represent fitting uncertainties, and 
do not include the error from the correlation with reddening. 



As shown in Figure 20, our best-fit solution yields (m — 
M) = 18.34 ±0.06 ±0.16 (P-L zero point). Extinction- 
corrected distance moduli at the best-fit reddening are listed 
in Table 18. The first error in distance is a quadrature sum of 
a fitting error (±0.05 mag) and an error from Ry (=f0.01 mag 
for ±0.5 in Ry). We also added in quadrature an error from the 
metallicity (=f0.03 mag for ±0.06 dex). As in the NGC 4258 
case the zero-point error of the Galactic P-L relations domi- 
nates the combined error in the LMC distance modulus. With- 
out metallicity corrections, we would derive (m-M)o =18.48. 

We also derived (E(B-V)} = 0.11 ±0.02, where the error 
is from the lcr contour in Figure 20. Our average redden- 
ing is 0.03 mag smaller than the average reddening from the 
OGLE-II, (E(B-V)) = 0.143, which was adopted in the HST 
Key Project. 11 However, our reddening estimate is closer 
to the value of E(B-V) ~ 0.10 (Walker 1999) from sev- 
eral different reddening estimates for LMC stars (see Ta- 
ble II in McNamara & Feltz 1980). It is also in agreement 
with the reddening value determined for the LMC Cepheids, 
(E(B-V)} = 0.076 (Caldwell & Coulson 1985; Laney & Cald- 
well 2007). The zero point of the OGLE-II reddening map 
was set by two LMC clusters (NGC 1850 and NGC 1835) 
and one eclipsing binary (HV 2274), but the reddening for 
HV 2274 has different values from various studies: OGLE-II 
adopted E(B-V) = 0.149 ±0.015 (Udalski et al. 1998), which 
is larger than E(B-V) = 0.120 ±0.009 (Guinan et al. 1998), 
0.088 ±0.023 (Nelson et al. 2000), and 0.103 ±0.007 (Groe- 
newegen & Salaris 2001). Furthermore, as shown in Table 18, 
extinction corrections from the OGLE-II reddening map re- 
sult in systematically shorter distances at shorter wavelengths, 
which may indicate an overcorrection for extinction (see also 
Fouque et al. 2003). 

Our distance modulus is 0.16 mag smaller (or 7% in dis- 
tance) than the adopted distance modulus by the HST Key 
Project (18.50 mag). Formally, this implies an increase in 
their Ho by 7% ± 8%. Our LMC distance is also shorter 
by 0.20 mag (or 9% in distance) than the distance modulus 
adopted by the Sandage et al. (2006) program (18.54 mag). 
In addition to our study, recent parallax studies from the HST 
(Benedict et al. 2007) and the revised Hipparcos catalog (van 
Leeuwen et al. 2007) used the Wesenheit index Wvi to derive 

11 The HST Key Project adopted the OGLE-II LMC P-L relations and 
scaled them to (m-M)o = 18.50. However, the P-L relations were not scaled 
to (E(B-V)) = 0.10 as they claimed (see also Fouque et al. 2003). 
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FIG. 20. — Left: Apparent LMC distance moduli in BVIcJHK s with the best-fitting CCM89 extinction curve (solid line). Middle: Likelihood contours in the 
average LMC reddening and its distance modulus in BVIcJHK, shown at A\ 2 = 2.30, 6.17, and 11.8 (68.3%, 95.4%, and 99.73% confidence levels for 2 degrees 
of freedom). The error bar represents the adopted distance and its error by HST Key Project at the average reddening from the OGLE-II reddening map. Right: 
Same as the left panel, but likelihood contours when BVI C photometry is used only. 



(m-M)o = 18.40 ±0.05 and 18.39 ±0.05 mag, respectively, 
also suggesting a downward revision of the LMC distance. 
Finally, Macri et al. (2006) estimated a distance modulus of 
NGC 4258 relative to the LMC of 10.88 ± 0.04 (statistical) ± 
0.05 (systematic). Given the maser distance, this was then 
translated into the LMC distance modulus of (m - M)q = 
18.41 ±0.10 (statistical) ± 0.13 (systematic), which is within 
1 a of our estimate. 

We note that stronger constraints on distance can be 
achieved when several photometric bands are used including 
near-infrared data. As illustrated in Figure 20, the statisti- 
cal error in distance becomes larger by a factor of 3 when 
only BVIc information is used (right panel) rather than from 
BVIcJHK s (middle panel). The error would be even larger 
when the distance and reddening values are estimated from 
only two bands, such as the Wesenheit index. 

6.3. M33 

As our final case, we estimated the average distance and 
reddening for M33 Cepheids using the same method as in the 
previous section. The distance to this nearby spiral galaxy 
has recently been measured from a detached eclipsing binary, 
(m-M)o = 24.92±0.12 (Bonanos et al. 2006), and from two 
water maser sources in H II regions, (m-M)o 
(Brunthaler et al. 2005). 

Figure 21 shows likelihood contours in the solution for the 
average reddening and distance modulus. As in the HST Key 
Project, we used BVIc photometry in Freedman et al. (1991) 
for their 10 best observed Cepheids, which was obtained at 
the 3.6 m Canada-France-Hawaii Telescope (CFHT), the Kitt 
Peak National Observatory (KPNO) 4 m telescope, and the 
Palomar 1.5 m telescope. We assumed 0.03 mag for the pho- 
tometric zero-point error in each band. On the Zaritsky et 
al. (1994) scale, M33 Cepheids have (12 + log(0/H)) = 8.82, 
which was adopted in the HST Key Project. We therefore 
applied metallicity corrections from Kochanek (1997) for 
A(log(0/H))(Galaxy-M33) = -0.02. 

At the best-fit reddening, we obtained an extinction- 
corrected distance modulus of 24.51 ±0.1 1, 24.60 ±0.09, and 
24.53 ± 0.08 in BVIc, respectively, where the uncertainties 
are fitting errors only. The average distance is (m-M)o = 
24.55 ± 0.28 ± 0.16 (P-L zero point). The first error was de- 
termined by adding in quadrature the fitting error (0.28 mag), 
the error from Ry (=F0-03 mag for ARy = ±0.5), and the error 
from the metallicity (q=0.08 mag for ±0.15 dex). We adopted 
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the error in the metallicity from the HST Key Project. If 
we adopt metal abundances derived from electron tempera- 
tures (Sakai et al. 2004), the abundance difference between 
the LMC and M33 becomes A(log(0/H)) = +0.21, which 
leads to A(log(O/H))(Galaxy-M33) = +0.09. We would de- 
rive (m-M)o = 24.53 without metallicity corrections. 

The HST Key Project used the same ground-based photom- 
etry for this galaxy and derived a metallicity-corrected dis- 
tance modulus of (m-M)o = 24.62 ±0.15 based only on VI 
as for the other HST sample galaxies. While this distance 
is in good agreement with our estimate, our average redden- 
ing of (E(B-V)) = 0.08 ±0.09 is lower than their estimate, 
(E(V-I)} =0.27 ±0.05 or (E(B-V)) = 0.21 ± 0.04. 12 A 
part of the difference is due to different zero points of P- 
L relations. However, our reddening is in good agreement 
with the value of the eclipsing binary, E(B-V) = 0.09 ±0.01, 
which was obtained from BVRJHK S photometry (Bonanos et 
al. 2006) and with previous estimates based on several differ- 
ent methods (see Table 7 of Bonanos et al.). 

A similar result was obtained from a larger Cepheid sample 
in Macri et al. (2001a). They provided BVIc photometry for 
25 1 Cepheids obtained at the Fred Lawrence Whipple Obser- 
vatory (FLWO) 1 .2 m telescope and the MDM 1.3 m telescope 
(see also Mochejska et al. 2001a,b). From the restricted sam- 
ple of 90 Cepheids in 1 .1 < logP < 1.7 similar to those of the 
Freedman et al. sample, we derived (m-M)o = 24.67 ± 0.18 
and (E(B-V)} = 0.04 ±0.06, where the uncertainties are fit- 
ting errors only. 

Distances from these two independent sets of photometry 
are about 0.3 mag shorter than the distance from the eclipsing 
binary, although our distance has a large error. One potential 
explanation for the short Cepheid-based distance is the un- 
resolved blends by a strong star-to-star correlation function. 
Mochejska et al. (2001c) investigated the influence of blend- 
ing by comparing HST WFPC2 and Macri et al. ground-based 
images for a sample of 102 Cepheids (see also Mochejska et 
al. 2000). They found that the average flux contribution from 
unresolved luminous companions on the ground-based im- 
ages could be on average <~ 20% in BVI for logP > 1, which 
would systematically underestimate the Cepheid distance by 
~ 10%. In addition, Lee et al. (2002) obtained single-epoch 

12 We note that Freedman et al. (1991) originally derived (E(B-V)) = 
0. 10 ± 0.09 and (m-M) = 24.64 ± 0.09 using BVRI photometry, assuming 
the LMC distance modulus of 18.5 mag and its mean reddening of 0.10 mag. 
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FIG. 21. — Left: Apparent distance moduli in BVIc with the best-fitting CCM89 extinction curve (solid line). Right: Likelihood contours in the average 
reddening and distance modulus of M33 Cepheids shown at A\ 2 = 2.30, 6.17, and 11.8 (68.3%, 95.4%, and 99.73% confidence levels for 2 degrees of freedom). 
Two open circles represent the distance and reddening estimates from an eclipsing binary (Bonanos et al. 2006) and those from the HST Key Project. The 
horizontal lines represent a distance from the water maser sources and its la error (Brunthaler et al. 2005). 



/-band photometry for a subset of the Macri et al. sample us- 
ing HST WFPC2 and found that the HST photometry is on 
average w 0.2 mag fainter than Macri et al. photometry. 

Follow-up observations with a higher spatial resolution 
(e.g., Macri 2004) would be helpful for a definitive test of 
the distance to M33, although there are indications that our 
estimated reddening is more accurate than the value in the 
HST Key Project. Note that the effect of stellar blending is 
less significant in our distance estimates for NGC 4258 and 
the LMC. The ground-based photometry with the full width 
at half-maximum (FWHM) of ~ 1.5" would correspond to 
<~ 0.3 pc at the distance of the LMC, compared to <~ 7 pc at 
the distance of M33. The ACSAVFC FWHM of 0.09"(Ford et 
al. 2003) corresponds to ~ 3 pc at the distance of NGC 4258. 

7. SUMMARY AND DISCUSSION 

We have continued our effort to improve the accuracy of 
isochrones and distances derived from MS fitting. We ex- 
tended the previous Hyades-based calibration to the upper 
MS by constructing empirical color- T e ff corrections to match 
the observed Pleiades MS at the cluster's accurately known 
distance. We applied these empirically calibrated sets of 
isochrones to Cepheid-bearing Galactic open clusters to de- 
rive distances, reddenings, and Ry at the spectroscopic metal 
abundances and obtained Galactic P-L relations based on nine 
Cepheids in seven clusters. Our distance modulus of individ- 
ual Cepheids has an accuracy of <~ 0.08 mag (or 4% in dis- 
tance), which is compatible to those of recent parallaxes from 
the HST (Benedict et al. 2007) and the revised Hipparcos cat- 
alog (van Leeuwen et al. 2007). 

Using these P-L relations, we derived a distance to 
NGC 4258 and found that our Cepheid distance is in excellent 
agreement with the maser distance, supporting our distance 
scale from the MS-fitting technique. From BVIcJHK, pho- 
tometry we derived the LMC distance modulus of (m-M)o = 
18.34±0.06±0.16 (P-L zero point) after applying metallic - 
ity corrections. This is also in close agreement with distance 
estimates from the recent HST and Hipparcos parallax stud- 
ies. However, our revised LMC distance is lower than the 
distance adopted by the HST Key Project, which formally im- 
plies an increase in their Ho by 7% ± 8%. A similar size of an 
increase in Hq is expected from the SN la calibration program 
by Sandage et al. (2006). 

Our reddening estimates are systematically lower than 
those in the HST Key Project. Part of the reason is the P-L re- 



lations with different assumptions about the reddening. They 
adopted the color excess for LMC Cepheids from the OGLE- 
II reddening map, but its average value is A(E(B-V)) ~ 0.03 
larger than our LMC reddening estimate. Nonetheless, rel- 
ative distances between the LMC and target galaxies in the 
HST Key Project would remain less affected because the 
reddening-free or Wesenheit index adopted in the HST Key 
Project is designed to avoid the problem of knowing the ab- 
solute zero point of the reddening scale. In other words, the 
absolute reddening value would not be important as long as 
the difference in reddening between target Cepheids and the 
calibrating system is well-defined. Therefore, the Wesenheit 
index can be effectively used to build the cosmic distance 
scale once highly accurate distances to calibrating Cepheids 
are available. 

However, our E(B-V) estimates for NGC 4258 and M33 
are larger by (E(B-V)} <~ 0.1 than those in the HST Key 
Project, more than expected from the difference in the LMC 
reddening. This is due to two reasons. First, the Cepheid's 
color is not only affected by reddening, but it can be also af- 
fected by the metal content. The HST Key Project corrected 
a distance modulus for a metallicity effect after deriving a 
color excess. On the other hand, we applied metallicity cor- 
rections to apparent distance moduli from Kochanek (1997), 
and then derived a reddening value. Because the sign and size 
of the metallicity corrections could not be the same in differ- 
ent broadband filters, the difference in the Cepheid metallicity 
could lead to different reddening estimates from those in the 
HST Key Project. Second, the Wesenheit magnitude is based 
on only two broadband filters. However, reddening and dis- 
tance estimates are naturally correlated in the standard pro- 
cedure to estimate these parameters using P-L relations, so a 
small zero-point error in the photometry, for example, can be 
translated into a large error in reddening. Therefore, obtain- 
ing Cepheid photometry at least in three or more passbands 
including near-infrared filters is of great importance in the 
distance and reddening estimation (e.g., Gieren et al. 2006; 
Soszyhski et al. 2006). 

In the determination of individual Cepheids' absolute mag- 
nitudes, the extinction is one of the largest sources of error 
because of strong differential reddening in Cepheid-bearing 
young open clusters. Although we found an error of ~ 
0.03 mag in the mean E(B-V) for a cluster, most Cepheid- 
bearing clusters have patchy reddening that makes the individ- 
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ual Cepheid reddening uncertain by up to <Je(b-v) ~ 0.12. An 
average reddening of cluster members located near a Cepheid 
was often used to determine the Cepheid's reddening, but dif- 
ferential reddening over a small scale could not be resolved 
unless we have a sufficient number of nearby stars with well- 
determined reddening. However, it is noted that the most se- 
cure reddening estimates for Cepheids are considered to be 
given by the cluster Cepheids. In addition, our LMC redden- 
ing is in close agreement with the previous standard value, 
(E(B-V)) ~ 0.10 (Walker 1999), and our average redden- 
ing for M33 Cepheids is in good agreement with the value 
inferred from BVRJHK S photometry of an eclipsing binary 
(Bonanos et al. 2006). 

In terms of the zero point of the Galactic P-L relations, the 
uncertainty in the metallicity scale of the Galactic open clus- 
ters is the largest error source, with a change of A[M/H] ~ 
0.1 dex producing a ~ 0.1 mag shift in distance modulus. No 
other source is unusually large. In Paper III we have shown 
that several color indices can be effectively used in MS fitting 
to determine a cluster metallicity because of differential sen- 
sitivities of these colors on metal abundance. Although the 
noise in the current photometry prevented the estimation of 
the photometric metal abundance, deep multi-color photome- 
try combined with near-infrared data will be useful to deter- 
mine R v , E(B-V), (m-M)o, and [M/H], simultaneously. A 
more accurate metallicity scale could then be used to reduce 
the size of the systematic error in the P-L zero point. In ad- 
dition, the zero-point error can be further reduced by deeper 



photometry of other Cepheid-bearing clusters, which were not 
included in this paper because of poor photometric quality. 
These are C1814-191, Collinder 394, Mon OB2, NGC 1647, 
NGC 6649, NGC 6664, Platais 1, Trumpler 35, Vel OB5, and 
Vul OBI. 
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